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Abstract 

We present the first complete calculation of the anomalous magnetic moment of 
the muon in the minimal Randall-Sundrum model with standard model fields in 
five-dimensional (5D) warped space and a brane-localized Higgs. We use a fully 5D 
framework to compute the one-loop matching coefficients of the effective theory at 
the electroweak scale. The extra contribution to the anomalous magnetic moment 
is 

Aa^ rj 8.8 • KT 11 x (1 TeV/T) 2 , 

where 1/T denotes the location of the TeV brane in conformal coordinates, and is 
related to the mass of the lowest gauge boson KK excitation by Mkk ~ 2.5 T. The 
result is robust against the variation of the bulk fermion masses and 5D Yukawa 
coupling. 



1 Introduction 



The idea [TJ that our four- dimensional world is one of the boundaries of a slice of strongly 
curved five-dimensional (5D) Anti-de-Sitter space provides an attractive approach to the 
gauge-gravity and flavour hierarchy problems of the Standard Model (SM). In conformal 
coordinates the metric of the 5D bulk is 



where k ~ Mj>\ ~ 10 19 GeV is of order of the Planck scale Mpi, while the metric on 
the four- dimensional boundaries located at z = 1/k and z = 1/T is flat. This set-up 
provides a solution to the gauge-gravity hierarchy problem, since, if the proper distance 
between the two branes, 1/k x \n(k/T), is only a few times the Planck length, then T 
can still be of order of the TeV scale. The original proposal [1] considered only gravity 
propagating in the bulk, but it was soon realized that the SM fields may be 5D, too [2H6]. 
An exception is the Higgs field, which in the simplest set-up is required to be localized 
on the TeV brane z = 1/T to solve the hierarchy problem. Allowing the SM fermions 
to propagate into the fifth dimension opens up the possibility to explain the quark and 
lepton mass, and the quark mixing-angle hierarchies [HEllT] without postulating large 
hierarchies or small coupling constants in the 5D theory. 

The model predicts Kaluza-Klein (KK) excitations of the SM particles with typical 
KK mass Mkk ~ 2.5 T for the first excitation. The collider physics and flavour phe- 
nomenology of these KK excitations has been extensively studied. Direct searches at 
the Large Hadron Collider lead to lower limits on Mkk similar to those for Drell-Yan 
like single production of heavy resonances, which are now in the 1 to 2 TeV region. In- 
direct constraints from electroweak precision observables and flavour require higher KK 
masses, but may be avoided by adding custodial protection or flavour structure. The 
overwhelming number of studies is based on investigating the implications of tree-level 
production or tree-level exchange of the KK excitations. 

Comparatively little is known about the Randall- Sundrum (RS) model at the loop 
level due to the intricacies of quantum field theory in curved space-times. It appears that 
contrary to 5D theories with a flat compact extra dimension, the running of gauge cou- 
plings is only logarithmic [SHEE]. Divergences associated with the fifth dimension show 
up in divergent KK sums. The interpretation of these divergences is not as straight- 
forward as in Minkowski space. A precise regularization scheme, required to be able to 
calculate finite renormalizations, has not yet been developed. At the same time, there are 
a number of processes of much interest for searching for deviations of the SM, which arise 
only at the loop level; among them are the anomalous magnetic moment of the muon, 
and electromagnetic dipole transitions in the lepton (/i — > cy) and the quark sectors 
(b — > sj), which have been considered in the RS model in [12] and [T31417] . respectively. 
Higgs production through gluon-gluon fusion has been discussed in [T81 - I20] . 

The systematics of renormalization and the calculation of loop processes is probably 
easier to approach from the 5D perspective P|[T5] instead of the KK decomposition of the 
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dependence on the fifth coordinate. In this paper we set up the propagator and vertex 
Feynman rules for the 5D quantum field theory in a slice of Anti-de-Sitter (AdS) space 
(in this part there is significant overlap with [15]) and then revisit the muon anomalous 
magnetic moment as our first application. There is only one previous calculation of 
(g — 2) M in the RS model [12], which is based on summing over towers of KK excitations. 
Since (g — 2)^ is not ultraviolet (UV) sensitive, it is well suited to demonstrate the 
workings of the 5D formalism without the complications of divergences from the Planck 
scale. The main result of the present paper is a complete calculation of (g — 2)^, since, 
as will be explained, our result differs considerably from |12j . 

The outline of the paper is as follows. In Sec. |2]we set up notations by defining the 
5D theory and explain our strategy for calculating the anomalous magnetic moment. We 
match the 5D theory to an effective 4D theory, in which the scale T is integrated out, 
keeping the relevant SU(2)xU(l) invariant dimension-6 operators that must be added 
to the SM. After expressing the Lagrangian in terms of the fields in the electroweak vac- 
uum, we calculate the non-standard contribution to (g — 2) M in the 4D effective theory. 
The matching coefficients of the dimension-6 operators are computed in Sec. |3J We work 
with 5D fields and avoid multiple sums over KK modes by integrating out directly the 
fifth dimension. We identify two relevant contributions, related to the electric dipole 
transition and a four-fermion operator. The actual integrations must be performed nu- 
merically. The result in the single-flavour approximation is shown in Sec. SJ We also 
discuss the dependence of the result on the bulk mass and 5D Yukawa coupling and the 
generalization to the case with three lepton flavours, and we compare our result to pre- 
vious ones. We conclude in Sec. [5j The Feynman rules of the 5D theory are summarized 
in Appendices \M and [Bj and in Appendix O we provide the explicit expressions of all 
diagrams for reference. 

2 From the 5D theory to (g — 2) M 

In this section we define the 5D theory, and then discuss the 4D effective Lagrangian at 
scales below T relevant to electroweak and flavour physics when Mew T. We then 
calculate the muon anomalous magnetic moment 

«, = ^ = 4 M + ^ ( 2 ) 

in the effective theory, that is, we express the non-standard contributions Aa M in terms 
of the matching coefficients of higher-dimension operators in the effective Lagrangian. 

2.1 The 5D theory 

We adopt the simplest set-up for the SM in the RS geometry ([I]), in which all fields except 
the Higgs doublet can propagate in the bulk, and no further structure is added. While 
not the most attractive model for phenomenology since the scale T must be quite large 
due to the lack of custodial protection, it provides the appropriate starting point for the 
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study of loop processes. Quarks and the entire strong interaction sector are not relevant 
to leptonic processes at one-loop. In the following we specify only the SU(2)xU(l) 
interactions and leave out the quark sector. 
The 5D theory is defined by the action 

s m = J d'x 1^ dzVcl- \f mn f mn - \w a > MN w MN 



i— B! r. n L 



+ 

tjj=E,L 



GF+ghost 



+ J d 4 x |(W)tD M $ - V(Q) - (j) [y^iL^Ej + h.c] |, (3) 

where G = l/(kz) 10 is the determinant of the 5D metric and e^f = fczdiag(l, 1, 1, 1, 1) 
the inverse vielbein0 The potential is given by 

v($)=-^ 5D) ^y$t$+^($t$) 2 . ( 4 ) 

The covariant derivative is defined as usual by 

D M = d M - ig' 5 jB M - ig 5 T a W a M . (5) 

Y is the hypercharge of the field and T a the SU(2) generator in the given field repre- 
sentation (r a /2 with r a the Pauli matrices for the doublets Lj, $, and for the singlets 
Ei). The hypercharge and SU(2) field strengths read 

Fnm = d N B M - d M B N W^ M = d N W M - d M W% + g 5 e abc W b N W M . (6) 

The Higgs doublet is localized at z = 1/T and described by a four- dimensional field. 
The powers of (T/k) in the Higgs action arise from making the kinetic term canonical 
when the ^-integration in the original 5D action is eliminated using the 5(z — l/T) factor. 
The 5D lepton fields in the lepton-Higgs Yukawa interaction are therefore evaluated at 
z = 1/T. Latin indices from the middle of the alphabet refer to lepton flavour. 

The dimensionful parameters g5,g' 5 , M^^y^j of the Lagrangian are all of order M P1 . 
It is conventional to introduce the dimensionless quantities = M^Jk for the bulk 
fermion masses, whose precise values determine the hierarchical masses of the known 
leptons. The 5D action is given in the flavour basis that renders the bulk mass matrix 
diagonal, which can always be arranged by a field rotation. It is not possible in general 
to simultaneously diagonalize the 5D Yukawa coupling matrix y^^- Unlike the SM, the 
5D theory violates lepton flavour conservation and CP-symmetry in the lepton sector. 

The Feynman rules for the Higgs propagator and vertices involving the Higgs field 
can be read off directly from ([3]) in the standard way. The rules for the 5D fields in the 
bulk of AdS are discussed and summarized in Appendix |A] 



lr The spin connection term in the fermion action does not contribute and has been dropped. The 
gauge-hxing action is discussed in Appendix [X] 
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2.2 Integrating out the KK scale 

We assume that the scale T is much larger than the masses of the SM electroweak gauge 
bosons and the Higgs boson. In this case we may integrate out the fifth dimension, that 
is, the towers of KK states to obtain an effective four- dimensional theory. Since the states 
with masses below T are precisely represented by the SM fields, the four-dimensional 
theory is the SM plus SU(3) xSU(2) xU(l) invariant higher- dimension operators. 

An issue that arises here is that the 5D theory is non-renormalizable and must itself 
be defined as an effective theory below a scale A that should be at least a few times 
the Planck scale. It is generally believed that in the mixed representation the four- 
dimensional loop momenta should be cut-off at a value A(z) that depends on the position 
z in the fifth dimension. If A(l/k) is a few times the Planck scale, then the cut-off A(l/T) 
relevant to processes dominated by physics near the TeV brane should be near the TeV 
scale. This would render the naive application of the 5D formalism problematic, since it 
encodes the sum over all KK states rather than including only the few below the cut-off 
A(l/T). However, for a finite quantity such as the anomalous magnetic moment of the 
muon, the KK sum must converge, and the effect of including the entire tower relative to 
the truncation is of order T 2 /A(l/T) 2 , which is the generic size of corrections expected 
from the UV completion of the RS model. 

The dominant non-standard effects in the expansion in M-ew/T are associated with 
dimension-6 operators. We write 

i 

Since the matching coefficients q are dominated by distances C 1/T, the Higgs bilinear 
term in V($) can be treated as a perturbation, and the q can be computed in the theory 
with unbroken electroweak gauge symmetry. 

The effective Lagrangian is a special case of the general SU(3) xSU(2) xU(l) invari- 
ant theory with SM field content including dim-6 operators |21j . Since the anomalous 
magnetic moment is generated only at the one-loop level, we distinguish two types of 
operators. 

Class-0 — operators in C e e that contribute to at tree-level, but whose short- 
distance coefficients are generated in the RS model only at the one-loop order. The only 
relevant terms are 

Ci °*\ «*»-o = a B ,ij LiQa^EjB^ + a w>ij L l r a ^a^ u E j W a ^ u + h.c. (8) 

i 

Note that we use the same notation for the 5D and 4D fermion and gauge fields. It is 
understood that the fields in the 4D effective Lagrangian correspond to the zero modes in 
the KK decomposition of the 5D fields. The general effective Lagrangian also contains the 
chirality changing operators (LD^E)D^& and (LD ^E)D fl &, but they do not generate 
anomalous magnetic couplings of the photon at tree level. After replacing $ by its 
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vacuum expectation value, the operators are (L 2 d fl E)[D tl }22 and (L 2 d fl E)[D' I }22- Since 
[-D M ]22 contains only the Z-boson, this contributes to the anomalous magnetic couplings 
of the Z-boson but not the photon. The dimensionless coefficients cib,w are of order 
g^Vij^k 5 ^ 2 , where here g 5 denotes a generic 5D gauge coupling^ 

Class- 1 — operators in £ e fr that contribute to a M only at the one- loop level, but which 
are generated in the RS model at tree-level. The relevant terms are 



E 



+ c 2il {La^i&iD^) + c 3ii {UYr-L^ir-D^) , (9) 

where ir a D^ = 1/2 (ir a D^ — %D ^t"). They are generated at tree-level by the exchange 
of the 4D vector components of the KK gauge bosons. Although the 5D gauge cou- 
plings to fermions are flavour-diagonal and universal, flavour-dependence of the coeffi- 
cient functions arises through the different zero-mode profiles of the 4D fermions. The 
operators (LiD^E^D^Q and {LiD fl Ej)D^& discussed already above might contribute 
to the photon anomalous magnetic moment at the one-loop level, but inspection shows 
that it cannot be generated at tree-level from integrating out the KK scale. Another 
class- 1 operator that could possibly contribute to g — 2 is L^E^^ (and its hermi- 
tian conjugate). It can only be generated at tree level by diagrams with three Yukawa 
couplings. In scenarios with an (approximately) brane localized Higgs these diagrams 
are non-vanishing only when one includes the "wrong-chirality" Higgs couplings, which 
come with an additional suppression factor of order M| w /T 2 [22]. We generally neglect 
such higher-order terms in the M^w/T expansion. KK gauge boson exchange also gen- 
erates further four-fermion operators with only left-handed and only right-handed fields, 
but these do not contribute to the anomalous magnetic moment. The mixed left-right 
lepton-quark operator (LE)(QU), on the other hand, does contribute at one loop, but 
it is not generated in the RS model. Other class- 1 operators that one might write down 
contribute to a M only indirectly, through a modification of the relation of the SM param- 
eters to those of the 5D theory. When expressing the SM contribution to a M in terms of 
SM parameters, this effect is already accounted for. Thus represents the complete 
list of relevant operators. The coefficient functions b, 01,2,3 are of order g\k. In particular, 
we note the presence of chirality-preserving fermion-Higgs interactions proportional to 
5D gauge rather than Yukawa couplings. 

We now express the effective Lagrangian in terms of the fields in the electroweak 
vacuum, in which the SU(2)xU(l) gauge symmetry is spontaneously broken to electro- 
magnetism. Thus, we replace 



2 The 5D gauge couplings are related to the standard 4D couplings by = ( \ n (k/T) ) 



1/2 



5 



and 



D ^ d ^ teQA ^ l ^T 3 -sl / Q)Z^- i ^(T 1 + tT 2 )W;-t-^=(T 1 -tT 2 )W^ (11) 

where cw (sw) is the cosine (sine) of the Weinberg angle, with c\y = g 2 /(g 2 + g' 2 )- 
For the fermions it is convenient to transform to the mass basis, in which the effective 
Yukawa interaction —i/ijLi^Ej of the 4D fields is diagonal. Tree-level matching results 
in the 4D Yukawa coupling 



^=g) 8 /l?(i/r)/^(i/r)yr 



1 ~ 2C ^ / 1 + 2C ^ ..(5D) fc {m 

^ 1 _ {T/kf- 2 ^ V 1 - (T/k) l+2cE i Vij 

where f^\z) are the fermion zero- mode profiles, see appendix El The transition to the 
mass basis is given by the unitary rotations 

Li — >■ UijLj, Ei — > VijEj (13) 

of the 4D fields in £ eff . To go to the broken theory, 

Ei VijPRiPj, Li U^Pl^ ) , (14) 

where ^ is the Dirac spinor field for the massive leptons (i = 1,2,3 corresponding to 
electron, muon, tau) and v\ is the corresponding neutrino spinor field. Pl/r = (lT7s)/2 
are the chiral projectors. Inserting these substitutions into (ED, ([HD, generates numerous 
different operators composed of SM fields. We neglect operators which cannot contribute 
to the magnetic and electric dipole form factors or only start contributing at two loops. 
The effective Lagrangian for the computation of the anomalous magnetic (and electric) 
moments, as well as charged lepton-flavour violating processes at one loop is then given 
by 

CiOi -> Re (ai j )-^='ipi<j flu 'ilj j F flu + Im (aij)—^= ipia^il^jF^ 

+ ii,a \ (^Pli^^h) + [ l2M + l3>ij ) v - {faPtaMiiWH) 

+ 73,y ^$iPiaM)(-idr<lr) + 7a« ^(^P^^OM^) 

+ h.c. of previous line (15) 
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with Ffj, v the photon field strength tensor. The couplings are 

aij = [U ] aV}ij, 

ftijkl ^ [U^i m U m j\V^~\k n V n l b mn , 
~ ^ [V^irnVinj 



7i 



rn 



lfx,ij — 2. j [U^]imU m j Cx,m — 2,3) (16) 

m 

with aij = cw<iB,ij — Swa>w,ij- Many further terms can in principle appear on the right- 
hand side of ( I15p upon replacing the fields in Oi by those of the broken theory. Inspection 
shows that these additional terms are not relevant, since they contribute to the form 
factor Fx only, see ( TTT|) below, or to F 2 and the anomalous magnetic moment only 
beyond one loop. 

2.3 The muon anomalous magnetic moment 

The muon-photon vertex function is given by 



r M (p,p) = ieQputf ',s') 



7^i(g 2 ) + ?^ 2 (g 2 ) + 



2m M 



u{p,s) (17) 



where the ellipses denote parity- violating terms, and q = p' — p. The factor ieQ^ with 
Qn — —1 corresponds to our definition of the gauge covariant derivative and the standard 
normalization of the charge form factor -Fi(O) = 1. The anomalous magnetic moment is 
defined by 

= = W (18) 



It is straightforward to compute the extra contribution Aa M from the effective Lagrangian 
( IT5|) . Since a is a one-loop quantity, while (3 and j x are tree-level, we need to compute 
the diagrams shown in figure [TJ 

The tree-level diagram (a) with an insertion of the class-0 electromagnetic dipole 
operator results in 

r*W) = y* Re (a 22 ) 2iu(p', s')ia^q v u{p, s). (19) 

The imaginary part of the couplings generates an electric dipole moment. 

There are two possible contractions of the four-fermion operator (diagram (b)), which 
differ by the exchange of Pl -h- Pr and (3 2 kk2 Pk22k- The loop integral is ultraviolet 
(UV) divergent by power counting and must be regularized. The divergence is related 
to factorization of the contributions from the KK scale T and the low-energy scales me, 



7 



(b) (cl) (c2) 

73 p> J) 73 p> p S ~\73 p' p 73/'' ^ 



(dl) (d2) (el) (e2) 

Figure 1: Diagrams for a M from (j!5p . 

which is the only scale present in diagram (b) from the external momenta and massive 
lepton propagators. We adopt dimensional regularization with d = 4 — Is. The final 
result is finite, and arises from the UV 1/e pole that multiplies a numerator of 0(e). 
It depends on the scheme and we use the "naive dimensional regularization" (NDR) 
scheme with anti-commuting 75. The scheme dependence of this ultraviolet-sensitive 
term is compensated by a corresponding dependence in the calculation of the matching 
coefficients as we discuss later. The situation is similar to the calculation of the 6—7-57 
transition with the weak effective Lagrangian method. 

The contribution to the vertex function from the four-fermion operator in the effective 
Lagrangian is 

i/3 2 kk2 . n 2s f ddl «(P'> s')yP L i(J + m/j7"z(/ - 4 + m ih )^ p P R u{p, s) 



(2nY (P-ml)((l- q y-ml) 

+ (/?2fcfc2 <-> I3k22k, Ph ++ PR) 

= -ieQ^ - u(p, s')ia^ ' q v [P 2 kk2 + Pk22k + {^2kk2 - Pk22k)l5] u{p, s) (20) 

where a sum over internal lepton flavours k — 1, 2, 3 is understood, and the limit d — > 4 
has been taken. The 75 term is related to the electric dipole moment; the remaining 
terms are the desired contributions to the anomalous magnetic moment of the muon. 

The Higgs derivative interactions from ( [15]) contribute via the diagrams (cl), (c2) 
and (dl), (d2) of figure [TJ Diagram (cl) with an insertion of the 72 coupling is given by 

r"(p,p') = %r-e,| 
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2£ f d d l (+t)u(p>,s')P R {t(f -l + m tl ) 1 Hty-l + m ll )u(p,s) 

(I 2 - m 2 H )((p - I) 2 - ml)((p' - I) 2 - ml) ' K ] 

Here = \f2m ll jv is the small standard model muon Yukawa coupling. The prefactor 
is proportional to the muon mass m M = y^v/y/2, since one of the vertices is the SM 
Higgs-fermion vertex. We may therefore neglect the muon mass in the integrand, and 
use the on-shell conditions to obtain 

1 {P ' P) ~ 2T 2 l6Q ^ J (2k)* (I 2 - m 2 H )(l 2 - 2p • I) ' W 

This vanishes, since the loop integral can only result in $u(p, s) = 0. Diagram (c2) and 
the insertions of the 71,3 interactions vanish in an analogous way. The reason for this 
can be deduced without calculations from the operator ^PL/ppf^d^H , which, up to a 
total derivative, can be converted into d^^PL/R^^H . In the NDR scheme, by the field 
equations, this operator is proportional to lepton masses, or small coupling constants, 
and hence of higher-order in our approximations. The corresponding arguments cannot 
be applied in the unbroken theory, since there the operator contains two Higgs fields. 

Diagrams (dl) and (d2) involve an internal charged Higgs and a neutrino line. For 
diagram (dl) we find 

r*(p, P = ^(*e)Hy,) 

2£ f d d l % 2 u(p', s') P R (jl - fn + m u )P R u(p, s) (21 -p- p'Y 
* J (2k)* ((p - I) 2 - ml)((p> - I) 2 - ml)(P - ml) 1 } 

The prefactor is proportional to the muon mass, hence we can neglect the muon and 
neutrino mass in the remainder of the expression. Introducing the Feynman parameter 
x (x — 1 — x) to combine denominators leads to 

rW) = 

2e f f 1 . u(p',s')P R u(p,s) (p + p' -2l)» 



(27r) d J ((I — xp — xp') 2 — m 2 ^ + xp 2 + xp' 2 — (xp + xp') 2 ) 2 

" J w Jo 



e73,22«V ,, 2£ f d d l f 1 dx u(p')P R u(p)[(x - x)p^ + (x- x)p'» 



T 2 r J (2n) d J (I 2 -m 2 + xxq 2 ) 2 

= 0, (24) 

where in the last step we used that the integrand is antisymmetric under the exchange 
x ■<->■ x. For later purposes, we note that the zero arises due to a cancellation between 
a finite term and an ultraviolet-sensitive term that arises as the product 1/e x e of a 
UV divergence and a numerator of 0(e). To see this, we assume I 3> p,p', expand 
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the integrand of f l23|) to first order in the external momenta, and extract the (p + p'y 
structure that is relevant to the magnetic moment. This results in 



d d l u(p', s')P R u(p, s) (p + pf - 2iy ( 2p ■ I + 2p' ■ I 



^73,22^ p 2e 



T 2 (Au) d / 2 
xu(p',s')P R u(p, s) 



1 4\ r(e) , /NU 4r(i + £ ). 



= t: ^^HP,s')P R (p + pru(p,s) x g x i-i) (25) 

The total contribution is, of course, zero; but the contribution due to the UV pole is 
finite despite the superficial divergence of the integral. This will be important when 
verifying the scheme independence in Sec. I3.3.3I An analogous result holds for (d2). 

Diagrams (el) and (e2) follow from the insertion of the operator tpiP^^Vj (eA^^) 
in ffl5|) and its hermitian conjugate, respectively. None of the two contributes to g — 2 
as the loop only depends on a single external four- momentum. We can use on-shell 
condition to convert each appearance of to a lepton mass. The only possible Lorentz 
structure is then u(p', s'y^Pji/iuip, s)A fl , which contributes only to the Fx form factor. 

The muon anomalous magnetic moment follows by adding the non- vanishing matrix 
elements of the dipole and four-fermion operator, ( TT9|) and (|20|) . respectively, resulting 
in 

V ^ 7 1 o q 



4m; / Re (a 22 ) . 1 m/ 



T 2 \ y..e ' 167r 2 m 

\ yM fe=l,2,3 

Here = 1, 2, 3 refers to electron, muon and tau leptons, respectively. In passing to the 
last line of ( l26|) . we used that b mn , being generated by gauge interactions, is real. It 
therefore follows from the definitions (1TB]) that (3k22k = &2kk2i an d hence (3 2 kk2 + ftkiik is 
real as it should be. 



3 Calculation of matching coefficients 

We now turn to the calculation of the matching coefficients in the RS model ([3]). The di- 
agrams relevant to the anomalous magnetic moment in the full 5D theory with unbroken 
gauge symmetry are shown in figure [2j 
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Bla Bib B2a B2b 




B5a B5b 




Figure 2: Diagrams contributing to Aa^, more precisely the matching coefficients of the 
class-0 operators. Solid lines refer to leptons, with the right external line belonging to 
the doublet Lj, the left one to Ej. Wavy lines denote hypercharge gauge bosons and the 
external photon, curly lines SU(2) W-bosons. A solid-wavy (solid-curly) line refers to the 
scalar fifth component of the gauge field. Dashed lines denote Higgs bosons, including 
the external Higgs field (grey box). Vertices involving Higgs fields are localized at 1/T, 
all other vertices are integrated over position in the fifth dimension. 
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In a direct calculation of a M in the 5D theory, the diagrams should be represented 
in terms of the fields and interactions in the expansion around the electroweak vacuum 
of the spontaneously broken theory, which depend on the scales k, T, M EW , m ik . In 
drawing the diagrams in figure [2] we already accounted for the fact that we only need to 
extract the contributions from the scales k and T, which gives directly the short- distance 
coefficients of the SU(3)xSU(2)xU(l) invariant effective Lagrangian. The diagrams in 
the unbroken theory do not represent the infrared physics from the scales near and below 
the electroweak scale Mew correctly, but the incorrect infrared contribution cancels in 
the matching procedure. The low-momentum contributions from diagrams with 5D 
non-zero-mode propagators have already been included through the calculations of the 
previous subsection. Note that we do not consider diagrams with graviton exchange; 
see [12] for a discussion of the graviton contribution. 

The diagrams are understood in the mixed representation with an integration over 
four- dimensional loop momentum I and the vertex positions in the fifth dimension. A 
generic diagram such as diagram Bla in figure [2] can have three different contributions: 

• All 5D propagators propagate the zero mode. The loop integral does not contain 
the short- distance scales T, k explicitly, and is purely long-distance. This corre- 
sponds to a contribution to the SM anomalous magnetic moment a^ M . There is no 
need to perform this computation here, so we subtract this contribution. 

• At least one of the 5D propagators propagates a KK mode, but the loop momentum 
I satisfies ! <T. In this case, the subgraph consisting of propagators with KK 
modes can be contracted to a point. The reduced loop diagram must be computed 
with the 4D effective Lagrangian. This corresponds to tree-level matching of class- 1 
operators and the one-loop operator matrix element calculations performed in the 
previous subsection. An explicit diagrammatic analysis shows that the relevant 
subgraphs in figure |2] correspond precisely to the operators and neglected 
operators of dimension higher than six. In diagram Bla, for example, the only 
relevant contribution occurs when the gauge boson is a KK mode and all internal 
fermions are zero modes, which corresponds to the insertion of a local four-fermion 
operator as in diagram (b) of figure [TJ The fermion-Higgs effective vertex and 
diagram (dl) are related exclusively to diagram W8. If the gauge boson and at 
least one of the fermion propagators refers to a KK mode and the loop momentum 
is small, the contribution is suppressed and corresponds to a higher-dimension 
operator. If the gauge boson is the zero mode, then the flatness of the gauge zero 
mode and the orthogonality condition for the fermion modes forces all fermion lines 
to propagate zero modes as well. This is the SM contribution discussed above. 

• The loop momentum is of order I ~ T or larger. This contribution goes into the 
one-loop matching coefficients of the class-0 operators. 

In the following we describe the computation of the matching coefficients mentioned in 
the second and third item. 
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Figure 3: Hypercharge boson exchange that generates the four-fermion operator. 



3.1 Four-fermion operator 

The tree-level matching calculation is very similar to the one that generates flavour- 
changing four-quark operators from KK gluon exchange [131 [23-26]. Here we perform 
the calculation in the 5D picture, which makes it particularly simple (provided the 5D 
propagators are known). 

The relevant diagram is shown in figure [31 Only hypercharge vector gauge boson 
exchange can generate the (L i r y tJ, L i )(Ej'y^Ej) operator at tree level, since SU(2) gauge 
fields do not interact with E^. The scalar fifth component of the 5D gauge boson cannot 
be exchanged, since the two external zero-mode fermions at the vertex have the same 
handedness. From the tree diagram shown in figure [31 we obtain 

\2, 



% = -» Yf r />* w W AiMS(? = °' *• (27) 



Here Yl = —1, Ye = —2 denote the hypercharges of the leptons, and Aj_ is the r/^ 
component of the 5D gauge boson propagator given in H136[) of the appendix (see also [9]). 
We subtracted the zero mode from the propagator, since this corresponds to a SM 
contribution to a M . 

Because of the large mass of the KK excitations, the external lepton momenta of 
order m M can be set to zero. The zero-momentum limit of the gauge-boson propagator 
(HH is 

A ± (q, x, y) ^° eOn - y) ^ ( - i + \ j 1 ^^ ~ ^ ~ V 2 + ^ 

+ 2y 2 ln(yT) + 2y 2 In ^ J + 0(q 2 )^j +(x<r>y), (28) 

which agrees with a similar expression obtained in [25] from the explicit summation over 
all KK excitations. The 1/q 2 term is evidently the zero mode. Subtracting it, neglecting 
1/k 2 relative to 1/T 2 in the first term in curly brackets, since e = T/k ~ 1CT 16 is very 
small, and employing the normalization flSBJ) of the zero modes, we obtain 

1 /2 Yl Y e 2 



n Y 2 2 \T 2 In | 



13 



r l/T f^ 2 ( r ) fl/T a ^ )2 (y) 

/ dx^A/x 2 (2\n(xT)-l) + / dy "-ff^- y* (2 ln(yT) 



+ 2 In i [ UT d X dy ^ min I . (29) 



The remaining integrals over the fifth-dimension coordinates are dominated by x, y ~ 
1/T up to terms suppressed by some power of T/k provided the bulk mass parameters 
satisfy cz,. < 3/2 and c E , > —3/2, which will be assumed. We can therefore set the lower 
integration limit to and use the explicit form of the fermion zero modes from (]90aj) . 
(]90bp in the appendix to find 

kj =b + &i(cx.) + bx(-c Ej ) + b 2 (c Lz ,c Ej ) (30) 

with 

b = — - 



8 ln(l/e)' 
7 / \ ff' 2 (5-2c)(l-2c) e 2 - 1 

&l(c) = ~T (3 -2c) 2 rr^r» ^ 

, , , 9' 2 (1 - 2c L )(l + 2c £ )(3 - c L + c E ) . 1 e 2 ^- 1 e' 2 ^ 1 

bl\C L , C E ) = — — ry- ; In ■ 



4 (3 - 2c L )(3 + 2c E )(2 - c L + c E ) el- e 2 ^- 1 1 - e' 2 ^ 1 ' 

Substituting this result for b^ into (1T6|) and exploiting the unitarity of the flavour rotation 
matrices to the mass basis, we obtain 

Pw = boSij5 k i + 5 kl ^ b\{c Lm ) P% m U mj + c% ^ bi{~c En ) \V ] ] kn V nl 

m n 

+ $> 2 (cL m ,c B J [U% m U mj [V% n V nl . (32) 



m,n 



The first term does not change lepton flavour number F, the next two may change 
AF = ±1, and only the term in the second line can produce AF = ±2 transitions. All 
terms are formally required to compute Tlk m ^kl m ^ x fakk2 relevant to the anomalous 
magnetic moment. 

A particularly simple result is obtained in the "single-flavour approximation" , where 
we ignore the presence of other leptons than the muon. This corresponds to replacing 

Y,k m £k/ m v x fokk2 -> /?2222 = P, and Uij = Vij = 5ij, such that 

/2 ^ 

P = b + bi(c L ) + bi(-c E ) + b 2 (c L , c E ) = — — - /(e, c L , c E ). (33) 

8 ln(l/ej 

Setting further cl = —c E = 1/2 + r (r > 0), we find 



e 2r 1 2r(2 - r) 
1 — e Zr e (1 — rY 



e 2r , 1 
In - 



1 -e 



2r 



2 4r 2 



;i - r)(l - 2r) 



(34) 
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Figure 4: Weak isospin boson exchange that generates the fermion-Higgs operator. 



According to (|12p . r is related to the muon mass by 

2re 2r /,„ V2m, 



e 2r y^kv 



(35) 



With e = 10 16 and k = 1, we find r = 0.0749 and / = 0.95. In general, when 
yff^k is of order 1 and e within a few orders of magnitude near 10 -16 , the deviation of / 
from 1 is only a few percent, and the flavour-conserving four-fermion matching coefficient 
is dominated by the r- independent term 6 - in this case the total contribution (I3"3"j) is 
suppressed by the large logarithm ln(l/e) and the small hypercharge gauge coupling. 

3.2 Fermion-Higgs operator 

For completeness we also give the matching coefficient c^i for the fermion-Higgs operator 
(Zj7 M r a Lj)($tiT a _D At $), which gives rise to diagram (dl) in figure [TJ Only SU(2) gauge- 
boson exchange shown in figure H] can generate an operator with a single charged scalar 
at tree level. The fifth component of the 5D gauge boson cannot contribute due to the 
boundary conditions. We find 



c 3>4 = -i %^T 2 dx !^Af" s { q = 0, x, l/T) . (36) 



Inserting ([28]) into fl36l) we obtain 



\2. 



^ = -Y\^^-T^j llk (37) 

As above, we neglect 1/k 2 relative to l/T 2 , and keep only the leading term in e. Using 
the normalization of the fermion modes simplifies the expression to 



_9 2 

c 3,i — "3" 



lnl/e 



;i-2c L J(5-2c L J 2(1 - 2c Li ) In 1/e 



(3 - 2c L J 2 (3 - 2c 



Li , 



e 2cL 



1-e 2 ^" 1 



(38) 



Note that c 3ii is about 10 2 times larger than 6^-, since it is proportional to the SU(2) gauge 
coupling and since its largest term does not feature the suppression factor j— j-^ ~ 1/35 
that is present in 6y. For cl { larger than but not too close to 1/2 the dependence of c^i 
on the model parameters is rather mild, since the 1 in the round brackets dominates. 
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3.3 Electromagnetic dipole operator 

The matching coefficient of the electromagnetic dipole operator is more difficult to com- 
pute, since it comes from the 5D one-loop diagrams in figure El We are interested only 
in the electromagnetic dipole transition, so we take the external gauge boson to be the 
superposition cwB^ + sw W 3 '^ corresponding to the photon. We also anticipate that the 
Higgs doublet $ in the operators L^a^EjB^ and LiT a $a liV Ej W a,fiu will be replaced 
by its vacuum expectation value and therefore extract only the weak isospin component 
that corresponds to charged lepton — > transitions after electroweak symmetry 
breaking. These simplifications have already been employed to omit some additional 
diagrams with the photon coupling to the Higgs-doublet field. To compute the matching 
coefficient we apply the following strategy: 

• Subtract the zero mode from every internal 5D gauge boson propagator. It can be 
shown that selecting the gauge-boson zero modes forces the fermion propagators 
to only propagate zero modes. This follows from the flatness of the gauge boson 
and orthogonality of the fermion mode functions. Thus, the subtracted terms 
correspond precisely to the Standard Model contribution to g — 2. 

• Expand the diagram in the external fermion momenta p, p' and neglect the (tachy- 
onic) mass in the SU(2) Higgs propagator^ This ensures that we pick up only the 
contribution from loop momenta I ~ T as appropriate for a matching coefficient. 
Some of the one-loop diagrams have tree-level subgraphs corresponding to propa- 
gating KK excitations while the loop momentum is of order of the electroweak or 
muon mass scale. These contributions are excluded by the above expansion; they 
have already been accounted for by the other dimension-6 operators with coefficient 
functions b, q and the one-loop diagrams in Sec. 12.31 as discussed above. 



3.3.1 Sample diagram 

We now illustrate this procedure by discussing the calculation of diagram Bla from 
figure [2] in further detail. With the 5D Feynman rules from appendix |A] we find for this 
diagram the expression 

-iT 3 \ {5D) Y L Y E f 1/T dx f 1/T dy f 1/T dz f dH 



Bia = 0<? 5 ) («.<W (—J y, yy j i/k ^ J i/k m J i/k w J ] 



L t (p')P Rlp A^(p' - I, z, y) 7 „Af (p - I, y, l/T)Af (p - 1, 1/T, x^PrE^p), (39) 

which contains the 5D propagators of the fermions and the zero-mode subtracted hyper- 
charge gauge boson, and the wave functions f^' \z)L i {p l ), (x)Ej(p), f^\y)e*^ of the 
external states, which are the zero modes of the 5D fields. The positions of the three 



3 Recall that the computation of the matching coefficients is done in the unbroken gauge theory. 
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vertices in the 5th dimension, x, y, z, and the four-dimensional loop momentum I are 
integrated over. The electromagnetic charge prefactor arises, because the coupling of the 
external photon in isospin space is [-^1 + = Q» = ~ 1- 

We decompose each fermion propagator into its four chiral components using (|12ip . 
Most of the 64 possible terms vanish due to the projectors Pr in and the brane 
boundary conditions g^il/T) = /^.(l/T) = 0, see (IBBI) . (1891) . The two remaining terms 
can be deduced from (I123al) . which results in 



Bla 



g ,2 e 5 Q,Y L Y E y£ D) T 3 [V T dx f^ T dy f^ T dz f dH 



4P J 1/k {hxf J 1/k {kyY Ji/k (kz)* J (2tt)4 

fS(z)fi%)g^(x)^A^(l,x t z) 

W) *> tf)*£(P. V> l /T)F^{v, 1/T, x) {l P (p"-/h,l„} (p - I) 2 + 

d + F^(p, z, y)d-F+(p, y, l/T)F^(p, 1/T, x) b P l^-^} ] P R Ej(p), (40) 

where p = p—l, p' = p'—l. In the KK picture the first term in square brackets corresponds 
to picking up three momentum factors from the propagator numerators, while the second 
contains two KK mass factors. In appendix Owe provide the expressions corresponding 
to fl39|) . (j40p for all 21 one-particle irreducible diagrams shown in figure [2j 

We note that at this point the integral over the coordinate y of the external photon 
vertex could be performed analytically. Since the photon zero mode is constant and 
combines with e^f\°\y) = e to the dimensionless 4D electric charge, the remaining y- 
integral is 



I" 1 dy v H)Wl!/!! 
Jx/k \ky) „ iP 2 -rnl){p 2 -ml) 



and a similar expression for the mass terms. The key point is that the ^/-integral is 
the orthogonality relation of the fermion mode functions. Thus, the KK number is not 
changed at the external photon vertex and the remaining sum involves only two mode 
functions. The KK sum can be expressed in terms of Bessel functions. However, the 
result is algebraically complicated, so we do not make use of this simplification in the 
numerical evaluation of the diagrams. 

Returning to (|40p we perform the expansion in the small external momenta p, p' and 
pick up the terms linear in the momenta, since these produce the electromagnetic dipole 
structure ia^ v q v after application of the Gordon identity. Since the various F-functions 
depend only on the square of the four- momentum argument, this expansion can be done 
by using, e.g., 

Fl{p' - /, z, y) = F+(l, z, y) - 2p> • I -^F+(l, z,y) + .... (42) 
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The derivatives of all required propagator functions are listed in appendix [B] In Feynman 
gauge (where the gauge boson propagator is proportional to rj pl '), the first term in square 
brackets in (j40p turns into 

F+(l, z, y)F+(l, y, l/T)F^(i, 1/T, x) {l v fl^ v I 2 + 7 77„7, 2p-l) 

-2p-lF+(l,z,y)^[F+(l,y,l/T)F^(l,l/T,x)]\, (43) 

where terms odd in the loop momentum I have already been dropped, since they integrate 
to zero. The angular integrations can now be done trivially. In the above integral we 
simply replace l a l@ — > r] a/3 l 2 /4: and perform the Dirac algebra. Employing the on-shell 
conditions, 

i't'lvlv 4 P ;, 777^7, 2p ■ I 2p M Z 2 (44) 

etc. What remains is a scalar integrand that depends on I 2 only, and the bulk coordinates 
x, y, z. The integral can be evaluated numerically after performing the Wick rotation 
I 2 — > —I 2 . The propagators given explicitly in the appendix refer to these Wick-rotated 
("Euclidean") propagators. Diagrams involving internal Higgs lines are easier to evalu- 
ate, since the number of independent bulk coordinates to be integrated reduces to two 
or even one, due to the brane localization of the Higgs field. 

For the numerical evaluation the loop integral needs to be ultraviolet and infrared 
finite. We find that the integration converges for large /, as expected since the leading 
one-loop expression for the anomalous magnetic moment should not require UV renor- 
malization. The diagrams are also IR finite with the exception of the internal insertion 
diagrams Bla and Bib, for which only the sum is finite, and the non-abelian diagram 
W8. The IR divergence arises when all internal fermion modes are zero modes. There 
is no IR divergence in the internal insertion diagrams B3a, B3b with scalar hypercharge 
boson exchange, since the chirality-flip at the fermion- B 5 vertex forbids the propagation 
of zero modes in these diagrams. The existence of an IR divergence in Bla, Bib and W8 
implies that the purely four-dimensional treatment described above potentially misses 
terms of the form 1/e x e. We discuss this further below. 



3.3.2 External insertions 

Another comment is necessary on the diagrams with an external Higgs insertion, such 
as B2a. The fermion propagator that connects the external Higgs vertex to the internal 
gauge vertex is (in this example) 

Af(p,x, 1/T)P R = -F+(p,x, l/T)jP R + cTF+(p,x, l/T)P R . (45) 

If this expression contained a relevant 1/p 2 contribution from the zero mode, the diagram 
would be long-distance sensitive, and the expansion in the small external momenta would 
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be invalid. The second term is a mass term, see f ll23a[) . to which the zero mode cannot 
contribute, since mo = 0. Hence this term can be expanded in p, and since it depends only 
on p 2 , we can simply set p = to linear order. The zero mode is present in the first term, 
which contains ip ! /p 2 x f^\x)f^\l/T). If the one-particle pole at p 2 = remains in the 
final answer, then this part of the external Higgs insertion into a zero mode needs to be 
amputated; it corresponds to the first term in the sum of tree diagrams that sums to the 
SM lepton mass matrix. After amputation the remaining short- distance contribution 
is a one-loop correction to the chirality preserving LiLtf vertex. The general L%La 
vertex function with off-shell zero-mode fermions can be decomposed into "on-shell" and 
"off-shell" terms as follows: 

A M = A ^ + /A^ p , + A ^ p ^. (46) 

The first piece constitutes the correction to the on-shell LiLa vertex and is not relevant 
to the anomalous magnetic moment. It constitutes the long-distance piece that must 
be amputated. The off-shell term with a f" vanishes, since diagram B2a has an on-shell 
external p' line on the right. However, the A^ ff p p* term cancels the propagator factor 
p'/p 2 from the internal fermion line, and represents a short- distance contribution that 
must be added to the matching coefficient of the dipole operator and hence contributes 
to g-2E 

It follows that in all external insertion diagram^ we have two different contributions 
to consider. The first one can be obtained by replacing 



Af(p,x, 1/T) dTF+(p = 0,x, 1/T) = —9 { S{x)ff{l/T) 



, fcV/ 2 {kx) l ' 2 - c ^ - (kxf^- 1 ' 2 

where the last expression follows from the explicit expression for the propagator functions 
given in the appendix. We note that the KK mode contribution to d~F^ (p, x, 1/T) is not 
suppressed due to the large KK masses of order T. Thus, the external insertion diagrams 
are not suppressed relative to the internal insertions, contrary to what has been assumed 
in the previous literature, where the external insertions have been neg lected^ 

The second contribution stems from the off-shell vertex terms and explicitly contains 
the zero-mode propagator. It can be obtained via the replacement 

A"Af (p,x, 1/T) -> A^J x % ijp 2 x f%\x)f%\l/T) = i A^ p ff(x)ff(l/T) , (48) 

where the internal zero-mode pole is cancelled. We will refer to these second contributions 
as "off-shell" terms. They require the separate computation of the A^ ff p structure of the 



4 See [37] for a related discussion in a different context. 

5 B2a, B4a, Wf, W2, W4, W5, W6a, W6b in figure H Diagrams B2b, B4b are treated similarly, 
with the appropriate modifications for an insertion into a line with momentum p' . 
6 See, however, the arXiv version [28] of [IB] . 
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LiLij vertex, which can be obtained from the expansion in the small external momenta. 
The expansion must now be performed to second order, since we need the ^p^, ^p ,fJl 
structures in the vertex function A M . 

The second contribution appears to be suppressed by a power of the lepton mass 
matrix due to the two additional fermion zero mode profiles. The expansion of A M to 
second order brings an additional factor 1/T from the scale of the loop momentum. The 
ratio of the second to the first contribution can therefore be estimated as 



TJ L 



(0), 



d-F+(p = 0,x,l/T) 



1 " 2c Li 



-i 2 



1 - (T/k) l - 2c ^ 



(49) 



where we put x = 1/T for the second estimate. For symmetric bulk mass parameters 
= —cl, such that /_£ (x) = g^(x), and anarchic Yukawa couplings this is indeed of 
order m^/v, see f[T2T) . For the other special choice cl — > 0.5 we find 



d-F£(p = 0,x,l/T) Txlog(kx) 



(50) 



For x close to 1/T this approaches — l/log(l/e) « —1/35, and we obtain a lepton-mass 
independent suppression. Thus, the "off-shell" terms are expected to be subleading 
for standard choices of the bulk mass parameters, which is indeed found in the exact 
numerical evaluation discussed in Sec. 14.11 



3.3.3 Scheme independence 

We mentioned above that the naive four- dimensional calculation of the sum of the indi- 
vidually IR divergent diagrams Bla+Blb might be wrong, since it could miss terms of 
the form 1/e x e. A similar issue appears in the finite result (1201) for diagram (b) from 
figure [H for which we would obtain zero, if the Dirac algebra were performed naively 
in four dimensions. The finite results of both calculations depend on the treatment of 
7 5 in d dimensions. In the following we show algebraically that the scheme dependence 
cancels if the same scheme is applied consistently to both parts, and that the finite result 
agrees with what we obtained above. 

Both parts are limits of a diagram in the full 5D theory with fields expanded around 
the electroweak vacuum and massive zero-mode fermions. This diagram is finite (as far 
as the a^ v q v structure is concerned). The divergences arise only when this diagram is 
split into a short-distance contribution from the KK scale, diagrams Bla+Blb, and a 
long-distance contribution from the electroweak and lepton mass scale, represented by 
the four-fermion operator insertion diagram. Before factorizing the full-theory diagram 
we can freely anti-commute 75. The convention that corresponds to (1201) amounts to 
eliminating all chiral projectors except for two placed at the operator vertex 7 p -Pl®7 p -Pr. 

For consistency, the same convention must be applied to the starting expression for 
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the diagrams Bla+Blb, which therefore reads 



I \2 



1T 3 



l/T 
1/k 



f(0)' 



k 3 J V2 

soy 



d d l 



i j K x ) 9e^ [ z j > - 

{kxY {kzY ^ J (27T 



i L m {p') lp P L — 



(p> - If 



Y 



7 



+ 



r 



(p - I) 2 (p' - I) 2 1 J (p - Z) 



-/ 



■ lv P R E n {p). (51) 



To obtain this result, we write down the explicit expressions from appendix [C] and 
multiply it by [t/t] 2i ^ 2 so that the external states correspond to the SM muons in the 
mass eigenbasis. All fermion propagators can be replaced by zero-mode propagators, 
since only this term is IR divergent, and the chiral projectors are placed as discussed 
above. Performing the integration over the bulk coordinate y of the photon vertex, we 
arrive at ( IBTj) . 

The IR divergence comes from I T, hence we can set I = in the zero-mode 
subtracted gauge-boson propagator. Moreover, only the r\ pv structure contributes to 
the IR divergence. After these simplifications, we can identify the expression (J27l) for 
bij. We perform the expansion in external momenta and simplify the Dirac algebra as 
much as possible. When the Dirac structure is multiplied by a 1/e pole we do not anti- 
commute 75 with 7^, but since the external momenta are four-dimensional we may use 



that a(W 



75 



"75 J 



W. The result of all this is that the IR sensitive contribution in the 



matching coefficient of the electromagnetic dipole operator is given by 



leQn — nrTnn u(p , s ) 
IQ71 T 



2 (p + p'YPR - ~ YP^ u qul p P R 



u(p,s). (52) 



The second term in brackets cannot be further reduced without assumptions on 75. 
Since YPl^^^vIpPr vanishes in four dimensions, the entire expression is finite, but 
scheme-dependent . 

In the one-loop matrix element calculation of the four-fermion operator that led 
to (1201) we assumed the naive dimensional regularization (NDR) scheme to obtain the 
result in the last line. We repeat the calculation in an arbitrary scheme, leaving out 
the (/?2A:fc2 ^ /3fc22fc, Pl ** Pr) term, which is related to the matching coefficient of 
the electromagnetic dipole operator with L and E exchanged. It is easy to see that 
the ultraviolet 1/e pole in the one- loop matrix element (|2"U|) that is relevant to the 
electromagnetic dipole transition is obtained correctly by approximating (J + m£ k )'~f fl (J — 
$ + rri£ k ) — > —mi k Yi- Dropping terms proportional to g M , we obtain 



r uv(P>P') = ie Q^ -^^ ^u{p',s') YPdo-^q^pPR u(p,s). 
Now we can make the following two observations: 



(53) 
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• In the NDR scheme ( ant i- commuting 75), between u(p', s')[. . ]u(p, s), 

YP L i^ v q»l P P R = 2e(p + p'YPr = -2e ia^q v P R , (54) 

where the Gordon identity is used in the last step, and the term not related to the 
electromagnetic dipole structure is dropped. Then ( 152]) is zero and fl53|) coincides 
with f l20|) . Thus we reproduce the previous results. 

• In an arbitrary scheme the scheme- dependent terms drop out when summing 
rf R (p,p') and T^j Y (p,p'), and the result agrees with the one obtained before in 
the NDR scheme. 

This proves that the result is scheme- independent, and that the naive four- dimensional 
treatment of diagrams Bla, Bib fortuitously gives the correct result. 

A similar issue arises for the IR-sensitive diagram W8, but in this case it turns out that 
the naive four- dimensional calculation needs to be corrected by a finite term of the form 
l/exe that is present when working consistently in 4 — 2e dimensions. This finite term is 
scheme-dependent in general, and the scheme- dependence is compensated by the effective 
theory diagram (dl) in figure [TJ The zero result f )23|) for diagram (dl) corresponds to 
a full-theory diagram in which the charged-Higgs vertex in f jT5|) is written as ^i-PfOV^j 
rather than ^fipPiyy We must apply this convention consistently to diagram W8 as 
well. The IR sensitive part of W8 can be obtained from (I163p by replacing the internal 
fermion by the zero mode. Then, rotating the external states to the mass eigenbasis by 
multiplying (I163P by the relevant part of W8 simplifies to 



r'W) 







r 3 " 








2 _ 




f d d l 


rl/T 

dx 

l/k 




J (2nY J 


(kxy 



" f 2 ' l— d J dx J -i^£^tf-l,x,l/T)f!< l >(l/T) 

Exploiting q v A^usiQ, x i V^) = A^ MS (g, x, 1/T, £) , which follows from (11311) . 

T 3 ?/ 5D) ?; 

y^f* (im${im = X>* u lk t kh vi 3 , (se) 

and e&fy (x) = e, this simplifies further to 

1 r H d l r 1/T f^ )2 (x) 

r*W) = 9!(ie)- mix J2[U%U i2 ^ dx^-^A^(p>-l,x,l/T,0 

xL{p )Pr (p- W -iyp PrE{p) ■ (57) 
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Since the difference to the naive four- dimensional calculation arises only from I <C T, 
we apply a cut-off f C A C T to the loop momentum. This allows us to set the gauge- 
boson propagator momentum to zero without generating a spurious UV divergence, 
and to identify the coefficient function c 3)i , since A™ s (0, x, l/T, f) = A5_ MS (0, x, l/T). 
Expansion to linear order in p, p 1 results in 

A 

r>W) = (-e) £ c 3 Au%u i2 ^ / I 

xL(p')P R E(p) . (58) 

The difference between the correct and naive treatment comes from the second term 
in square brackets, where we must use IH V — > 1/d x Z 2 ?]^ instead of 1/4 as done in 
the four-dimensional treatment of the integrandQ Using the definition (TT6|) of 73,^, we 
therefore find that the difference between the d- dimensional and four- dimensional result 
that must be added to the naive four-dimensional treatment of W8 is 

A 

= ieQ» |g§| x I^PhCp + pO M ^(p) + 0{e) . (59) 

The result is independent of the arbitrary cut-off as it should be and is part of aw,ij } the 
coefficient of the class-0 operator L^^a^Ej W a '^ u n We can check f|59l) by noting that 
it should be minus the UV sensitive e xl/e term in the matrix element calculation that 
led to (T25T) . Since = — 1 this is indeed the case. 

In principle the issue of scheme-dependence could also occur for diagrams B5a, B5b 
and W3, related to the effective theory diagrams (cl), (c2). However, it can be checked 
that in this case our convention for the position of the chiral projection operators in ffT5|) 
implies that diagrams (cl) and (c2) are always zero and accordingly there are no e x l/e 
terms in B5a, B5b and W3. 



(p + p) 



4l^(p + p') ■ I 

T 2 



3.3.4 Gauge invariance 

To check the gauge independence of the matching coefficient, we performed the calcu- 
lation in 5D i?£ gauge, and verified that our numerical result is independent of £ for a 
range of values of £. We also checked analytically that the ^-dependent terms cancel. 
The proof is too lengthy to be presented here explicitly, but we find it instructive to 
discuss the structure of the argument and the key algebraic identities. 

7 In this case the square bracket vanishes, which explains the absence of an explicit IR divergence. 

8 One could choose a scheme where W8 agrees with the purely four-dimensional treatment. In this 
scheme we would find a non-zero contribution from the one- loop diagrams (dl) (and (d2)). Our choice 
corresponds to the NDR scheme. 
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It is clear that the subsets of diagrams with hypercharge and SU(2) gauge boson 
exchange must be separately gauge-independent. It is also evident that the five diagrams 
Bxa and the other five hypercharge exchange diagrams Bxb must be separately invariant, 
since the dependence on the bulk mass parameters c^. and Ce is different whether the 
Higgs coupling to the fermion is to the left or to the right of the one of the external 
photon. Hence we consider first Bla to B5a, see figure |2j 

The gauge-parameter dependence arises from the scalar gauge boson propagator A 5 
(see (I134p . d 138[) ) and the longitudinal part 

^ A,, (I, z, z', = ^ A±{l/Vt, z, z') (60) 

of the vector gauge boson propagator. In diagrams Bla, B2a and B5a we therefore 
replace A pu (l, z, z') by the previous expression. Using (11061) . the mode representation 
of the propagator, and the completeness relation ^ n fj%{z)ffy(z') = kz5(z — z') we 
derive the identity 

A B (l, z, z') = ^5(z- z') - i d z d z ,A x (l/^, z, z'), (61) 

which is used to eliminate the propagator of the scalar component of the gauge boson 
in diagrams B3a and B4a. At this point all gauge-dependence is in the expression 
Aj_(Z/-\/£, z, z'), which appears in all five diagrams. The delta-function term in ( IBTT) 
results in a ^-independent expression that can be ignored]^ 

To proceed we need an identity that relates the different diagram topologies - internal 
insertions (Bla, B3a), external insertions (B2a, B4a), and Higgs (B5a) diagrams. This is 
obtained as follows. After using ( )6T|) . we integrate by parts the bulk coordinate deriva- 
tives on Aj_, which yields the corresponding derivatives on the two fermion propagators 
(B3a), or one propagator and one external zero mode (B4a) adjacent to the gauge boson 
vertex. Now we can use successively identities such as 

d z (jj^u A i(Ph x ' zhsAf(p 2 , z, y)^j = Af (pi, x, z) - jfe) Af(p 2 , z, y) 

+ zAf ( Pl , x, z)5(z -y)- tAf(p 2 , z, y)5(x - z) (62) 

which follow from the defining equations for the fermion propagator fllOSp . (11091) . If one of 
the propagators on the left-hand side of (|62|) is an external zero mode, the corresponding 
delta-function and momentum term is absent. When this is used in B3a and B4a, the 
first term on the right-hand side of (|62l cancels precisely the gauge-dependent terms 
of Bla and B2a, respectively, leaving over only the delta-function terms and the Higgs 
diagram B5a. 

Some of the delta-function terms vanish, since the remaining integral is independent 
of p' (p), such that the loop integral must result in $E(p)j = (Li(p')p'' = 0). The 

9 This is not true in the non-abelian diagrams with more than one internal gauge boson line. 
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surviving two terms correspond to diagrams with one deleted propagator. Inspection 
shows that these two terms and the Higgs diagram B5a are identical up to a different 
hypercharge prefactor. The prefactors are such that we obtain 

Y L [B3a] - Y E [B4a] - [B5a] = -1 + 2-1 = 0, (63) 

which proves the gauge-parameter independence of the abelian diagrams Bla to B5a. 
The algebra works out in the same way for the symmetric diagrams Bib to B5b. 

The proof of gauge invariance of the SU(2) gauge boson contributions though similar 
is more involved and we only make a few remarks. First, the three abelian diagram 
topologies Wl to W3 are by themselves gauge-independent, and the proof proceeds as 
above. The remaining eight genuinely non-abelian diagrams W4 to W10 can also be 
split into two groups using the structure of the vertices and propagators. It is useful to 
begin with diagram W5 containing two internal scalar W bosons, with 7W5W5 vertex 
ip — l)n + (p' — I) ft, an d to proceed by treating the two terms separately. To show the 
gauge cancellation, identities such as 



yd y ~d y Aj_(k } y,x) 



f 1/T dy f 1/T dy 

/ TT [3Al(*> z >I/>0] d y A x (k,y,x) = - / — A\ { (k, z, y, f) 

Jl/k K V Jl/k K V 

(64) 

must be used, where the derivative term on the right-hand-side can subsequently be 
eliminated by the equation of motion fl!32a[) for the gauge-boson propagator, and the 
boundary terms in the partial integration vanish due to the orbifold boundary conditions 
on the branes. To simplify the proof of gauge-invariance of the genuinely non-abelian 
diagrams we assumed that the external photon has physical transverse polarization, 
e* ■ q = 0. 



3.3.5 Anapole moment and one-particle reducible diagrams 

The calculation of the diagrams of figure [2] does not produce the expected form factor 

(p + p'ru(p',s')P R u(p,s) } (65) 

which can be converted into the anomalous magnetic moment structure via the Gordon 
identity. Instead, we find different coefficients of p M and that is, there is an extra q 11 
term. More precisely, we find a nonzero coefficient C q for the structure L q^P^E arising 
only from the genuinely non-abelian diagrams W4 to W10 in figure [2j It can be shown 
that the same diagrams, but with an incoming doublet and outgoing SU(2) singlet field 
lead to the appearance of the structure E q^P L L with coefficient — C q . Thus, the vertex 
function ( TT7|) contains the parity-violating term 

r"(p,j/) D ieQ^C q u(p\s')q^u(p,s). (66) 

We find analytically that the overall dependence of C q on the external lepton mass is 

C q = 2m t Cf (67) 
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Figure 5: One-particle reducible diagrams with short-distance contributions. 



where C red is independent of the mass of external states. Therefore (I66p can be rewritten 



as 



r»(p,p') D teQ,Cl ed u(p',s')^ l5 u(P,s) 



(68) 



The presence of such a term does not necessarily violate the Ward identity for the vertex 
function. The most general vertex function of on-shell fermions compatible with U(l) em 
gauge invariance is given by (see, e.g., [29|[30] ) 



r"(p,p') =ieQ ll u(p',s' 



rF 1 (q 2 ) + 



io^ v q v 
2m, , 



F 2 (q 2 ) + 



2m,, 



l 5 F 3 (q 2 ) 



+ (gV-^)75^4(g 2 



u(p, s) 



(69) 



The form factor F4 is related to the so-called anapole moment of the lepton. Contrary to 
the magnetic and electric dipole moment the anapole moment is gauge- dependent [29] 
and by itself not a physical observable. It is, however, non-zero in the Standard Model 
[3Tj . Since we dropped the g 2 7 At 7s terms in the course of our calculation, the left-over 
q^ terms fl68|) should be interpreted as the RS contribution to the anapole moment of 
the muon. They are irrelevant for the dipole moment calculation. The (^75 term can be 
removed either analytically or by averaging the coefficients of p^ and p ,fl to extract the 
coefficient of {p + p'Y- 

In fact, there is an entire class of one-particle reducible (1PR) diagrams shown in 
figure O which generates only terms oc 5^75. The diagrams in the first row give a short- 
distance contribution if the internal Higgs propagator connecting to the fermion line is 
cancelled. In this case the contributions are finite and have exactly the same structure 
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as the term arising from the diagrams of figure [2] discussed above. The diagrams in 
the second row of figure [5] are different in that they are not only gauge dependent but 
also suffer from UV divergences. This is not a problem per se as the anapole moment 
is known to have both features already in the SM [21]. In any case, we do not have to 
consider the contributions from the self-energy like 1PR diagrams further, since we are 
not interested in the anapole moment. 



3.3.6 Numerical evaluation 

The evaluation of the Wilson coefficient of the electromagnetic dipole operator requires 
the numerical evaluation of up to four- dimensional integrals. The integration variables 
are typically three 5D coordinates and the modulus of the loop momentum. For I 3> T 
the Bessel functions in the propagators behave like exponentials and it is possible to 
determine the asymptotic behaviour of the integrand analytically. The numerical inte- 
gration of the integrand in the high-momentum region is time-consuming and introduces 
significant uncertainties. The origin of this problem is large cancellations between the 
different Bessel functions in the propagators. Thus, we perform the 5D position integrals 
over the whole (length) of the fifth dimension, but limit the integral over the modulus 
of the momentum to the region below some momentum cut-off. In the high-momentum 
region we replace the integrand by its leading asymptotic behaviour, and integrate it 
from the cut-off to I = oo. In practice a cut-off of O(100T) is already sufficiently high. 
In order to estimate the uncertainty associated with this approximation we vary the 
cut-off in the interval from 60T to 150T; the typical uncertainties are of the order of a 
few per mille. 

Some of the propagator functions feature discontinuous behaviour if two bulk coor- 
dinates coincide (see appendix IA.3j) . Both speed and precision of the evaluation can be 



improved by explicitly splitting the integration domain into sectors where the integrand 
is continuous. This approach appears to be slightly more efficient than using a coor- 
dinate transformation to make the discontinuities manifest by mapping them onto the 
coordinate axes. Separating the different regions leads to a reduction of uncertainties by 
about a factor of 2 to 10, depending on the amount of discontinuous propagators (and 
thus regions) in the diagram. 

The calculation is performed in general covariant gauge. Since the gauge invariant 
subsets of Feynman diagrams are known (see Sec. I3.3.4p a numerical verification of gauge 
invariance provides a powerful check for our code. Furthermore, the error estimates 
provided by standard integration routines can be checked by using the small spurious 
residual dependence of Aa^ on the gauge parameter £. Very large gauge parameters are 
numerically difficult to handle, but the limit £ — > oo, i.e. unitary gauge, can be performed 
analytically on the level of propagators and gives consistent results!^ For the error 
estimates we typically choose £ = 1 and £ = 4 as reference gauge parameters. This choice 



10 It is not generally true that the £ — > oo limit can be taken before integration over loop momentum. 
However, for the anomalous magnetic moment terms we find that they are already gauge-parameter 
independent after integrating over bulk positions, at fixed value of 4D loop momentum I. 
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already changes the contributions of individual diagrams significantly and provides a 
reliable check on the gauge-independence of the sum within numerical uncertainties. For 
instance, the result for the gauge invariant set of diagrams with an internal hypercharge 
boson (Bl-B5b) varies by one to three percent when changing the the gauge parameter 
by a factor four, while the individual diagrams experience changes of order one. A 
similar behaviour is found for the genuinely non-abelian diagrams (W4-W10); here the 
variation of the sum does not exceed two percent. The abelian W diagrams also form 
a gauge invariant subset. However, the sum of Wl, W2 and W3 is over two orders 
of magnitude smaller than the individual diagrams. This large cancellation makes the 
result numerically unstable and prohibits a check of the £ independence. Due to this 
cancellation the abelian W diagrams are not important for the determination of Aa M . 

It should be noted that the gauge invariance proof in Sec. l3.3.41 does not discriminate 
between "on-shell" and "off-shell" terms as it takes into account the full propagator of 
the external fermion. That is, only the sum of both sets of terms needs to be gauge 
invariant. Indeed we find numerically that the "off-shell" terms alone are ^-dependent, 
but we checked analytically that the gauge- dependent terms are the form (p + p')^ and 
hence can combine with the other terms to a gauge-independent result. Due to the 
relative smallness of the "off-shell" terms (see also Sec. 14.11) it is not possible to verify 
within the numerical accuracy that the on-shell terms contain a tiny residual gauge- 
dependence that cancels this ^-dependence. The analytic proof, however, shows that 
this must be the case. 

As a further check we work with two separate implementations of our evaluation 
strategy, which both rely on Mathematica to handle the numerical integrations. After 
combining the uncertainty due to the extrapolation and the error estimate for the in- 
tegration itself, we determine Aa M with a typical uncertainty at the percent level. The 
overall evaluation time generally depends on the specific choice of the 5D input param- 
eters and the desired precision; a typical runtime for the results presented in Sec. 0] is 
18 hours on an Intel i7-950 3.4 GHz processor. 

3.4 Comment on genuine Higgs contributions 

We note that the diagrams in figure [2] contain internal Higgs lines, but they do not 
correspond to the usual Higgs exchange contributions. For example, the SM Higgs 
contribution to g — 2 is not part of the zero mode contributions of these diagrams. 

Drawing a Higgs exchange diagram with an internal insertion, see figure E] left, we see 
that such diagrams cannot exist at the level of dimension-six operators L^a flu EF' J ' L/ , since 
the internal Higgs propagator would correspond to the (non-existing) contraction of 
the Higgs doublet, while the external line would be $^ rather than $. The external Higgs 
insertion diagram does have the required field content, but it contains chiral components 
of the brane-to-brane fermion propagator that vanish because of /£.(1/T) = g^il/T) = 
0. 

The leading Higgs exchange contributions are generated at one loop by diagrams such 
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Figure 6: Higgs boson exchange diagrams. The diagram on the left does not exist, 
since there is no propagator. The leading contributions come from higher-dimension 
operators generated by diagrams such as the right one. 

as in figure [6] right. This corresponds to dimension-8 operators of the form 

_L x [^y^y^y i^EjF^® , (to) 

and implies an additional v 2 /T 2 suppression relative to the gauge-boson contribution. 
The loop integral is infrared divergent due to the two (massless) Higgs propagators, 
which results in an additional logarithm of T, similar to the In (m#/m M ) in the SM 
Higgs contribution. 

In the present paper, we do not address dimension-8 operators. The operators gen- 
erated by Higgs exchange, despite being suppressed, may, however, be of interest for 
flavour-changing processes, since they depend on the flavour structure yy^y of Yukawa 
couplings, whereas the dimension-6 operators are proportional to a single Yukawa ma- 
trix ng. 

In this context, we also note that we do not take into account wrong-chirality Higgs 
couplings which arise when the singular brane-localized Higgs theory is defined carefully 
by a limiting procedure [22]. Because these couplings are themselves of order (v/T) 2 , 
these Higgs contributions should be considered together with the dimension-8 operators 
in the effective SM. 



4 Muon anomalous magnetic moment 
4.1 Single- flavour approximation 

We now discuss our result for the muon anomalous magnetic moment. We will mostly 
restrict ourselves to the single-flavour approximation, and discuss briefly below why it 
should be a good approximation to neglect lepton-flavour changing contributions. 

In this case the matching coefficients a = «22 (dipole operator) and (3 = ^2222 (four- 
fermion operator) defined in ( !T5|) are functions of k, T, the two bulk-mass parameters 
cl, ce-i and the dimensionless 5D Yukawa coupling = y^^k. We fix the Planck scale 
quantity k to 2.44 x 10 18 GeV [TJ, and determine ce from the value of the muon mass 

= 105.7 MeV through f[T2"j) . for given T, cl, and Y^. The other parameters that enter 
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the analysis are the Weinberg angle = sin 2 6w = 0.231, the electromagnetic coupling 
a em (M z ) = 1/128.94, and the Higgs vacuum expectation value v = 246 GeV. 

The contribution to the anomalous magnetic moment from the four-fermion operator 
follows from (I2BD and (1531): 

™J «em "ij 1 



[Aa M ]« = -j£ ^ = g^- ^ j^r <*) - 1.2 ■ 10-" x (1 TeV/T) . (71) 



In the last expression we adopted T = 1 TeV and symmetric bulk mass parameters 
cl = —ce = 0.578..., which leads to f — 0.95, and made the dominant dependence 
on T explicit. Thus, the four-fermion operator makes a negligible contribution to the 
anomalous magnetic moment compared to the dipole operator discussed below, mainly 
due to the l/ln(l/e) ~ 1/35 suppression. We find that / ~ 1 is always a good approx- 
imation, independent of the values of T and Y^, when the bulk mass parameters are 
near the symmetric value. This can be different when ci is larger than about 0.65, or 
smaller than 0.5, but (except for extreme choices of cl where ce becomes positive) / 
is never large enough to compensate the l/ln(l/e) suppression. Therefore we drop the 
four-fermion operator contribution from the further discussion. 

The contribution from the electromagnetic dipole operator is determined by five- 
dimensional loop integrals. For their evaluation we have to rely on numerical integration 
and must add (|5"9"|) to take the additional term from W8 into account. Before discussing 
the result, we provide an order-of-magnitude estimate of this contribution. Since the loop 
is dominated by momenta of order of the KK scale T, one might expect the dimensionless 
Wilson coefficient a 22 in (|T5|) to be of order e x a em /(47r). According to (|26|) this would 
result in Aa M ~ a em / (4?r) x (m^v)/T 2 1.6 - 10 -8 x (1 TeV/T) 2 , which is very large. This 
estimate is, however, too naive as it does not take into account the wave functions of the 
external 5D states. The photon zero-mode profile fllOOp provides a factor of 1/ ^\n(k/T). 
The conversion factors for the gauge couplings g$ — > g^p amount to ln 3,/2 (A;/T), since 
the one-loop diagrams involve three gauge couplings. Thus, we obtain an overall factor 
In (1/e). Furthermore, the external fermion zero-modes provide a factor of g^\x)f^\z). 
Since the integrals over the bulk coordinates are dominated by x, z ~ 1/T, we expect that 
a 2 2 is roughly proportional to g^(l/T)/j°' ) (l/T). Together with the vacuum expectation 
value v/y/2 in ( 126]) from the Higgs insertion, this counts as a factor of order m ,,. The 
contribution to Aa^ due to electromagnetic dipole operators is then estimated by 11 ! 

2 

Aa, » In (1/e) x Loop(c L ,c E ,y {5D \T,e) , (72) 

where the remaining loop factor should be 0(1). This estimate indicates that the con- 
tribution from the dipole operators is enhanced by a factor of the order In 2 (1/e) ~ 
10 3 compared to the contribution from the four-fermion operator; this would lead to 
Aa M ~ 2.4 • 10~ 10 for T = 1 TeV, which is only a factor of three smaller than the present 
uncertainty of a M . 



11 Note that one factor of m M arises from the loop diagram, the other from the definition of the F2 
form factor as coefficient of a^ v q v /2m^ rather than a^ v q v . 
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Figure 7: Left panel: Dependence of the different contributions to Aa M on the KK scale 
T: — [Aa M ]s (dark grey/red, dot-dashed), [Aajjy (light-grey /orange, dashed) and the 
sum (black/blue, solid). Right panel: Aa M (cj,, Y^/Aa^CL = —cr = 0.578425, F M = 1) 
for cl = 0.51, Y^ = 1 (dashed) and cl = — cr, = 10 (solid) as a function of T. 

For our numerical study it is convenient to split Aa M in contributions from an ex- 
change of hypercharge B bosons and W bosons 



and to study them separately. The left panel in figure [7] shows the T dependence of 
— [Aa M ] s (dot-dashed) and [Aa^]^ (dashed) as well as the sum (solid) for the parameter 
choice Cl = —ce and Y^ = 1. Note that the contribution due to hypercharge bosons 
B has opposite sign than the non-abelian contribution. The dependence of Aa^ on T 
follows almost precisely the 1/T 2 x ln(fc/T) scaling, see (172 p . Fitting the numerical T 
dependence shown in figure [7]to 1/T X x ln(fc/T) one finds a best fit for x = 1.997±0.04. 
The contribution of (|59|) to [Aa^]^ is typically at the level of 2.5%. This is consistent 
with the absence of a ln(fc/T) enhancement in (|38|) . The left panel shows the result for 
Aa^ for different choices of cl and Y^ normalized to Aa At (c^ = — ce,Y^ = 1). We see 
that the ratio is quite independent of the KK scale: the T dependence is practically 
universal. Also different choices of cl and Y^ change Aa^ only mildly. 

From figure [7] we conclude not only that the correction to g — 2 is quite insensitive 
to the specific choice for ?/ 5Z) ) and the mass parameters cl and ce, but also that the 
estimate f J72|) agrees fairly well with the numerical results if not for some amount of 
cancellation between [Aa M ]w and [AaJ^. Explicit numerical values for a representative 
set of parameters are shown in table [TJ where a reference value of T = 1 TeV was used. 

The robustness of the result with respect to the 5D parameters can be understood 
in the following way. Consider the contribution from terms that do not have a fermion 
zero-mode propagator. As already mentioned above each external (SM) fermion comes 
with the corresponding zero-mode profile; in our case and f L ^ and we expect these 
mode factors together with the Yukawa to combine roughly into a factor of the lepton 
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Table 1: Dependence of the contributions to Aa^ on mass parameter asymmetry and 
Yukawa coupling. T has been fixed to 1 TeV. The uncertainties include the numerical 
error from the integration, and the uncertainty from the extrapolation and residual £ 
dependence. 



mass. The dependence on the 5D parameters will then arise from the zero-mode sub- 
tracted fermion propagators. Since the KK modes are quite insensitive to the 5D mass 
parameters the dependence is, as observed, mild. 

Contributions from terms containing an explicit fermion zero-mode propagator, such 
as the so-called "off-shell" contributions or the IR sensitive part of W8 have two addi- 
tional zero-mode factors; for instance 



yj^imfr (i/T) . /f oo/f ( Z ) ~ ^ • /f ( Z ). (74) 

The product f^\x)f^\z) is very sensitive to the choice of the 5D mass parameters and 
one expects a much more pronounced dependence on 5D masses and Yukawa couplings. 
However, as already argued in Sec. 13.3.21 these terms will typically be suppressed and, 
hence, this effect is not seen in the total. 

We will illustrate this using the non-abelian diagrams as an example. Consider the 
two special fixed doublet 5D mass such that cl = 0.5 and the symmetric choice 

cl = —ce- In all non-abelian diagrams contributing to the La^ v EW^ v structure the 
actual loop is independent of eg; thus after fixing cl = 0.5 the remaining dependence 
on ce from the external lines and on Y# will always combine to a factor of the lepton 
mass rri£ - irrespective whether a zero- mode propagator is involved or not. For the 
symmetric choice cl = —ce the identity f^ix) = (x) suggests that the explicit factor 
fi ( x )fL ( z ) for zero- mode propagators should count as an additional power of rri£ if the 
loop is dominated by the region very close to the IR brane. In this case the term with a 
zero-mode propagator will scale like m| rather than as for the remaining contributions. 
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Figure 8: Lepton mass dependence of total non-abelian (solid) and non-abelian "off- 
shell" contributions to a M ; the dashed line refers to fixed cl = 0.5001, the dot-dashed 
one to the symmetric choice cl = —cr. Diamonds correspond to electron, muon and tau 
mass, respectively. 

In particular, the "off-shell" diagrams which by definition have an explicit zero-mode 
propagator are sensitive to the 5D mass parameters and not necessarily proportional to 
the lepton mass. 

To verify whether this intuitive picture is justified we study the lepton mass depen- 
dence of [Aa M ]iy numerically. Figure M shows the full contribution of all W diag rams 
compared to the contribution of "off-shell" terms as a function of the lepton masso Both 
T = 1 TeV and Yf — 1 are kept fixed. The total non-abelian contributions (solid lines) 
for both, the asymmetric and symmetric case are on top of each other and cannot be dis- 
tinguished on the logarithmic scale; both are proportional to mj. The "off-shell" terms 
show exactly the same scaling for fixed doublet modes with cl = 1/2 (dot-dashed line) 
and are a factor 60 smaller than the total contribution in agreement with our previous es- 
timate. For cl = —ce (dashed) the lepton-mass dependence is more involved. Following 
our argument above we expect mf behaviour. This is, however, only realized for lepton 
masses above 2 GeV. For leptons lighter than the muon we find same mf dependence 
as in the cl = 1/2 case. The source of this unexpected behaviour can be traced to two 
diagrams: Wl and W4. Both Wl and W4 contain a contribution where only fermion 
zero-modes propagate in the loop. Our argument for the expected suppression of the 
off-shell terms by one additional power of the lepton mass relies on IR-brane dominance 
of the coordinate integrals due to the localization of the KK modes (see Sec. I3.3.2p . For 

12 Only the sum over all diagrams is gauge independent; the pure off-shell terms are gauge dependent. 
For figures H] and [9] Feynman gauge was used. Furthermore, recall from the previous footnote that Aa M 
has one additional factor of me relative to the coefficient function aii- 
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Figure 9: Dependence of [a^w on cl- Squares correspond to the sum over all diagrams, 
circles to the "off-shell" terms only. 

the pure fermion zero-mode contribution this dominance is not present and these terms 
are only suppressed by a mass-independent factor instead of mi/ 'v. The zero- mode terms 
are thus less suppressed for very light leptons and give rise to the mj behaviour in the 
low mass region of figure [HJ The cl dependence of the different contributions to [a^] w is 
illustrated in figure IHJ Again sum over all diagrams shows almost no sensitivity whereas 
the small off-shell terms are quite sensitive to the specific choice for c^. 

Using the overall dependence on T 2 and the relative insensitivity to mass parameters 
and Yukawa couplings we can condense our result for the correction to g — 2 in the 
minimal RS model into 

Aa„ « 8.8 • KT 11 x (1 TeV/T) 2 , (75) 

which is our main result. The exact results differ from this approximation by less than 
10% for E (0.1,10), T e (0.5TeV,50TeV) and c L between 0.5 and the symmet- 
ric choice. The difference between the current experimental average of the anomalous 
magnetic moments of /z + and and the Standard Model prediction is given by [32] 

- af 1 = 287(63) (49) x 10~ n (76) 

where the two errors correspond to the combined experimental and theoretical uncer- 
tainties, respectively. Since we find Aa M to be positive in the RS model, a^ xp — ajp will 
be smaller than a^ xp — a^ M , decreasing the discrepancy between theory and experiment. 
However, the estimate f!75|) shows that for T > 1 TeV the correction to a^ M is too small. 
Even for a low KK scale of T m 500 GeV which corresponds to lowest KK excitations 
presently being excluded by direct searches, Aa M only reaches the level of 40 x 10~ n - 



34 



of order of the error on a M , but too small to reconcile theory and experiment at the 2a 
level. 



4.2 Beyond single- flavour 

The generalization to the case with three lepton flavours is in principle straightforward: 
the Yukawa factors have to be promoted to the 3x3 Yukawa matrix and the external 
muon states have to be represented in the mass basis, see ffTB"]) . In practice, this com- 
plicates the calculation, since it now depends on six bulk mass parameters and the full 
5D Yukawa matrix, and only three of these parameters can be eliminated by the known 
lepton masses. We now argue that this complication is not relevant for the anomalous 
magnetic moment. 

In the interaction basis the only source of flavour change is the Yukawa matrix itself. 
Thus, independent of its topology each Feynman diagram can schematically be written 

as 

J2ulMpy^M^fv jmi (77) 

%3 

where n = m = 2 for external muon states. Here M^ 1 ' and collect all factors depend- 
ing on only flavour indices i and j, respectively. For the single-flavour approximation 
we observed that the dominant contribution to the dipole-operator Wilson coefficient 
ctij (and thus to g — 2) is largely insensitive to the specific choice for the 5D masses or 
Yukawa couplings but depends linearly on the lepton mass. For the three-flavour case 
this means that the dominant contribution to g — 2 is roughly proportional to the mass 
matrix: 



M (1 V D M {2) 



«/ff(l/T)^W(l/r). (78) 

i 

As the matrices U and V diagonalize the mass matrix 



2^ Ulf^(l/T)yi; u> g^(l/T)V im = diag {m e , m„ m T } nm , (79) 



i,3 



the flavour dependence of the diagrams approximately combines into a lepton mass 
factor, irrespective of the specific choices of the bulk mass parameters cj_, v ce and 5D 
Yukawa matrix. Hence we can adopt the simplest choice for the calculation by taking 
the Yukawa matrix diagonal, which corresponds to the single-flavour approximation. 
The anomalous magnetic moment of the muon is therefore reliably determined in this 
approximation, up to small corrections to f!78|) . which are no larger than, e.g., the effect 
of unaccounted higher dimensional operators. 

The fact that (1771) is, up to small corrections, already diagonal, implies that lepton 
flavour-violating transitions are suppressed. Flavour-changing processes - such as /x — > 
e7 in the lepton sector or b — > 37 in the quark sector - are therefore sensitive to the 
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sub-dominant contributions that depend strongly on the 5D mass parameters. Such 
terms arise, for example, when the internal fermion lines of a diagram propagate zero 
modes and the gauge boson is a KK mode (see [H]), or via the contribution of the 
"off-shell" part of the LiL^ vertex with a subsequent mass insertion as discussed in 
the previous section. In this situation the suppression of "off-shell" diagrams is lifted 
and the terms have to be taken into account. Also Higgs exchange, though related to 
higher-dimensional operators, is relevant to flavour violation, since its different Yukawa 
matrix dependence produces terms that are not aligned with the mass matrix [T4~l[T5] . 

4.3 Comparison to previous work 

The first calculation of the muon anomalous magnetic moment in the RS model was done 
very soon [12] after the invention of the model. The authors used the 4D formalism and 
computed multiple KK sums. They (incorrectly) concluded that external mass insertion 
diagrams are suppressed by a factor (m^/M K K) 2 , which eliminates all diagrams except 
Bla and Bib in figure |2J In these diagrams a misplaced chiral projector at the photon 
vertex eliminates one of the two terms in the integrand f!148j) , f!149j) . The final result for 
the gauge-boson contribution to Aa^ in [12] is opposite in sign to ours, which is consistent 
with the fact that the abelian hypercharge contribution [Aa^]^ is negative. The final 
result f )75p is, however, dominated by M^-exchange, which was neglected in [12]. Ref. [T2J 
also considered the Higgs contribution and found that it puts important constraints on 
the bulk mass parameters. However, their Higgs contribution does not seem to contain 
the (v/T) 2 suppression relative to the gauge-boson contribution that is expected on 
general grounds, see Sec. 13.41 

The only other estimate of Aa M in the RS model with SM model fields in the bulk 
we are aware of is contained in [33]. The calculation is restricted to a one- loop vertex 
diagram with KK photon exchange but only zero-mode fermions in the loop. This subset 
of contributions significantly underestimates the true result, since it is suppressed relative 
to (172]) by a factor ln 2 (l/e) ~ 10 3 due to a cancellation in the photon KK sum. 

The calculations of lepton-flavour violating ^ — > processes [THUS] involve the 
same diagram topologies as g — 2, but the restriction to flavour violation allows to make 
some simplifying assumptions. Ref. [14J works in the KK picture and includes only 
the first KK excitation, while [15] adopts the 5D formalism as we have done in this 
paper. These works do not consider external insertion diagrams, though this has been 
corrected in [28]. The Higgs exchange diagrams that can be neglected for g — 2 provide 
the dominant source of flavour violation and are found to be of order 1/T 4 consistent 
with ( 170]) . Ref. [U] also considers the gauge-boson contributions, but assumes that only 
internal zero-mode fermions are relevant to flavour violation. Since we have seen that 
the KK fermions in the loop are by far the dominant contribution to flavour-conserving 
quantities like g — 2, it would be interesting (but numerically challenging) to investigate 
how the small non-alignment of these terms feeds into flavour violation. 
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5 Conclusion 



In this paper we performed the first complete calculation of the muon anomalous mag- 
netic moment in the minimal Randall- Sundrum model. We find that the additional con- 
tribution to a M is enhanced by a factor ln(fc/T) relative to the naive estimate a em /(47r) x 
(rn^/T) 2 , resulting in 

Aa M « 8.8 • 1CT 11 x (1 TeV/T) 2 , (80) 

nearly independent of the bulk mass and Yukawa coupling parameters of the model. For 
T « 500 GeV, which corresponds to first KK excitations around 1.3 TeV, this is of order 
of the theoretical uncertainty on a M . The minimal RS model is, however, already excluded 
by electroweak precision tests unless T is much larger, which makes the correction to 
g — 2 unobservable. T may be lower in models with custodial symmetry, but in this case 
the contribution from the extra states remains to be computed. 

The 5D framework developed here can be extended to an investigation of lepton and 
quark flavour- violating radiative processes in warped geometries that goes beyond the 
approximations adopted in [T4T4T6] . We hope to report on such results in the future. On 
the more formal side, it would be interesting to investigate the proper gauge-invariant 
regularization and renormalization of the RS model as an effective field theory, and the 
5D formalism seems to be the suitable framework. 
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A 5D Feynman rules 

In the following we provide the Feynman rules for the five-dimensional electroweak Stan- 
dard model ([3]) in the Anti-de-Sitter space time (pQ). The rules consist of the bulk mode 
functions, the vertices and the 5D propagators. 



A.l Equations of motion and mode profiles 
A. 1.1 Fermions 

The fermionic part of the action was given in ([3]). Upon inserting the explicit expressions 
for the metric and the vielbein we obtain for a generic fermion field ip 

1 1 4 



d x I dz 



kz 



+ iT\d z --)-- 

z z 



4 
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(81) 



where Y is the hypercharge of the fermion ip, and r a — > in case of a SU(2) singlet 
fermion. The dimensionless parameter c is related to the 5D mass M via c = M/k. The 
fifth Dirac matrix is chosen such that 



-175 



1 
-1 



(82) 



in the chiral representation. Splitting the fermion field into left- and right-handed com- 
ponents we obtain 



d x I dz 



4 



kz 



2 c - c - 

-^(02 )'0L ^i?^L V^i? ) + 

z z z 



where ipL = Pl^ and ^i? = Pr^P with 



R/L 



1±75 



(83) 



(84) 



and only the bilinear terms are given explicitly. 

Substituting the factorization ansatz ifj L / R (x,z) = f^ R (z)ip^ R (x) into the ac- 
tion, we obtain the mode equations [4] 



-/ R n) + ( 0: - ) /;; 



C („) 



(n) 



f (n) 



The orthonormality conditions 



/£°(*)/r(*) 



r l/T 


' 1 " 


/ 


ll/k 







' 1 " 


/ 


/cz 


h/k 



An) 



fR n \z)fP(z) = S 



(85a) 
(85b) 

(86) 



fix the normalization of the modes. The solutions for the zero mass eigenvalue are 
particularly simple: 



f R °\z)~z^, f?\z)~z*- c . (87) 

However, to obtain the SM in the low-energy limit we have to prescribe boundary 
conditions^ to exclude one of the zero mode of definite handedness. If ip describes a 
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This is equivalent to assigning different parities under orbifolding to left- and right-handed fields. 
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lepton SU(2) singlet, i.e. the right-handed massless field E in the SM, one needs to 
impose the boundary condition 

&\z) =0 (88) 

L V ' z=l/k,l/T V ' 

which removes the left-handed solution for the zero-mode and fixes the boundary condi- 
tion for the right-handed modes via (|85|) . Similarly, for left-handed SM fields the correct 
massless mode for the doublet lepton field L is obtained by imposing the condition 



= (89) 

z=l/k,l/T 



on the right-handed modes. 

In the following the labels on the fermion mode functions refer to the doublet (L) 
and singlet (E) field solution; in particular, the subscript L on a fermion mode or field 
will henceforth always refer to the lepton doublet L and not to its handedness. The left- 
handed (right-handed) mode functions will instead be denoted by f^(z) (g^(z)). The 
explicit expressions for the (massless) zero-modes are 



/i 0) (*) = y i-V/^ ^ (fcz)2(rg) " CL • (90b) 

The mode functions for non-zero eigenvalue are given by 

An), s n 5/2 / xr , \ Y c E +l/2 (^) Jc E +l/2{m n z) \ 

f E ( z ) = C E z' Y CE+1/2 {m n z) - , (91a) 

\ Jc E +l/2 { — ) J 

in), \ (~i 5/2 ( v I \ Yc E+l/2 (^t) Jc E -l/2(jn n z) \ 

9e (z) = C E z 0/ Y CE _ 1/2 {m n z) (91b) 

\ Jc E +l/2 { — ) J 



and 



An), \ fi 5/2 | \/~ , \ Y cl-1/2 (^) Jc L +l/2(m n z) 



(91c) 



in), s n 5/2 /v- t \ ^l- 1 / 2 (~T~ ) Jc L -l/l{m n z) . 

^ ; (^) = Cl^ 7 ( F Ci _i /2 (m n z) ^-v ) . (91d) 

Jc L -l/2 { — ) 



Applying the boundary conditions (188)) . (189)) . adjusted to the present notation, the mass 
eigenvalues m n are determined from 
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for the doublet and from 



Y 



/m r , 



,- L -\ 2 [—) 1 <•■>■; 2 



T 



-J, 



CE + l/2 



mr, 
k 



Y 



CE + l/2 



T 







(93) 



for the singlet field. The normalization factors Ce and Cl can be computed via the 
orthogonality relation (1861) . 



A. 1.2 Gauge bosons 

Let us start by repeating the expression for the different pieces of the action. The gauge 
part of the action for the non-abelian field can be written as 



Sr. 



d 4 x / dz 



1 



2kz 



W„ ( ri^d 2 - d»d v - rf v zd z ( - d z 



W u + 



- W 5 d 2 W 5 + 2W 5 d z d»W fl - 2g 5 e abc d^W^ h W^ c 

_ ?[ e abc e ade w b w c w »,d W »,e + g^^W^W^ 



(94) 



where d$ = d z , d 2 = and r/ 55 = —1 were already used and /x, v = 0, 1, 2, 3. The 
action for the abelian field B follows readily via the replacement W — >■ B and e abc — > 0. 

The gauge fixing part of the action is chosen such that the mixed W^W^ term in 
the second line of fl94l) is removed and no mixing of the 4D vector boson and the 
Goldstone W 5 can occur [9]: 



S 



GF 



d A x l dz 



-1 

1£kz 



d„W» - izd z I - W< 



The bilinear part of the gauge action then takes the form 



S 



G+GF 



d x { dz 



2kz 



1 



Wa d 2 ^ u - l - - d»d v - zd z -d z if v W„+ 



-W 5 d 2 W 5 + £W 5 d z [zd, 



Wf 



+ 



(95) 



(96) 



where the boundary terms left after integration by parts vanish as the W$ field vanishes on 
the branes. In the orbifold picture this follows directly if W$ is assigned odd parity under 
orbifolding; this is required to reproduce the observed low-energy particle spectrum. 
Covariant gauge fixing implies the existence of ghosts which, however, are not needed in 
this work, since ghosts do not couple to the W$ field. 
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Using the factorization ansatz W fl (x,z) = W^ n ^{x)f^\z) for the gauge field, 
the profiles fw(z) can be determined via the equation of motion 

zd z (~dj$(z)\=-m 2 n f$(z), (97) 



the boundary conditions on the branes 

An 
' W 

and the normalization 



d z f { w\z) 



= 0, (98) 

z=l/T,l/k 



^f { w\z)f { w\z)=*nm. (99) 
1/k Kz 

Note that the choice of the boundary condition follows from the need for a massless 
gauge boson mode to be present after integrating out the fifth dimension. 

The solution for zero mass eigenvalue has no dependence on the 5D coordinate and 
is given by 

Cm = (ioo) 

The mode functions for non-zero mass eigenvalue are given by 

f$\z) = C w z (jMnz) - ^f^ Y 1 (m n z)) , (101) 
where the eigenvalues m„ are given by the n-th zero of 

^- ^y - <io2) 

and Cw is determined by (I9"9"j) . Similarly, the mode equation for the Goldstone modes 
W5 is given by 

id z (zd a ±f$(z)) = . (103) 

The boundary conditions 

tiH= wt =° (104 > 

remove the zero-eigenvalue solution from the spectrum. One finds 

= C W5 z (j (m n z) - ^r^ Y (m n z)) (105) 

and the mass eigenvalues m n are the same as those for the 4D vector boson W^. 
Using the equations of motion yields the useful identities: 

d z f$\z) = m n f$(z) (106) 
dj$(z) = ^ 5 (z)-m n f%\z). (107) 
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A. 2 Vertices 



The Feynman rules for the vertices can be read off directly from the action. We denote 
non-abelian gauge fields W^ 1 by curly lines and the abelian field by wavy lines. The 
Goldstones, W5 and B 5 , are indicated by a superimposed solid line. 

Fermion-gauge field vertices 

In the unbroken theory all gauge interactions are flavour diagonal. Therefore, the flavour 
index % is conserved in all boson-fermion vertices. Below, a denotes the SU(2) index of 
the gauge boson of weak isospin and Y is the hypercharge for the fermion. For every 
vertex, there is an integral dz over the bulk coordinate of the vertex. 




Gauge boson self-interactions 

Here we only give the rules for trilinear gauge boson self-interactions. A, v and fi are 
Lorentz indices and a, b, c are SU(2) indices. The coordinate space derivatives always 
act on on the 5-coordinate of the vertex and only on the propagator /external state with 
the indicated SU(2) index. For every vertex, there is an integral Jy£ dz over the bulk 
coordinate of the vertex. 




.l_ g5e abc^ p _ q)v%x+{q 



l) ll V\u + {l-p)x r lu^ 
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o 
o 
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+ 1-95^ 

KZ 



aba 



jit, a 



v. b 



-ij-95 e abc Vvv {d z 

KZ 







z I on 'b' 



Higgs interactions 

As the Higgs field is IR brane localized all interactions occur at 1/T in the 5th dimension; 
therefore, we omit the 5-coordinate of the vertices. The only flavour changing interaction 
of the theory (with our basis choice) is introduced by the Yukawa coupling of the Higgs 
to fermions. Below m and n are SU(2) indices of the Higgs doublet, i and j are flavour 
indices. The Higgs field is represented by dashed lines. 
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wvwvw 



1 ' 2 \S1 11 

-g 5 Y rj^l 



. Yr a 

W5g 5 -rrv^ 



a I b 



43 



A. 3 Propagators 

The 5D propagators for bosons and fermions have been derived in the literature, see [9] 
and p3] respectively. However, we use slightly different conventions, e.g. our choice 
for the sign of the fermion mass term in the 5D Lagrangian differs from [15]. For 
completeness, we will therefore briefly sketch the strategy used to derive the propagators. 
Note that we try to use the same notation as the existing literature in order to allow for 
an easier comparison with previous results. 



A. 3.1 Fermion propagators 

The fermion propagator is given by the inverse of the 5D Dirac operator 

1 4 

VA(p,z,z') = i5(z-z')t. 



kz 



(108) 



In the mixed 4D momentum/bulk coordinate space representation we can replace i$ — > ft 
and obtain 



V = ft + iT\d z --)--. 

1 z z 



Following [15] we define a new function F via 



A(p, z, z) 



-ft-iT\d z 



z z 



F(p,z,z') . 



(109) 



(110) 



F is also a matrix in Dirac space and once we know F we also know the full fermion 
propagator. Inserting (IllOp into f ll08D and using the explicit representation for iT 5 gives 



-p 2 - d 2 + 



2 i c — c- 



-p 2 - d 2 + 



o 

,2 , c 2 +c-6 , 1 



F(p,z,z') = i(kz) 4 5(z - z')t 



4x4 



(111) 



where each entry of the matrix is itself a 2 x 2 matrix. The structure of (11 111) suggests 
a solution of the form 



F(p, z, z') = P L F-(p, z, z') + P R F + (p, z, z') 



(112) 



where F are now scalar-valued functions. The differential equation then can be solved 
in two separate regions: z < z' and z > z' . In each region one finds two equations 



-p 2 - d 2 + 



-p 2 - d 2 z + 



c 2 - c - 6 4 



+ -d z 

z 



c 2 + c-Q 4 



~d z 

z 



F-(p,z,z') = 
F + {p,z,z') = 



(113a) 
(113b) 
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which can be solved: 

F ± (p,z,z') = 



z 5 / 2 (A>J c±1/2 (pz) + B>Y c±l/2 {pz)) 
z 5 / 2 {A<J c±1/2 (pz) + B<Y c±1/2 (pz)) 



Here, J and Y are Bessel functions and the coefficients A± K and B±^ are functions of 



for z > z' 
for z' > z 

>>/< 



(114) 



z' and p that are fixed by the boundary conditions. Continuity at z = 
F±{p, z' + e, z') - F± (p, z' - e, z') -> for e^O 
and the "cusp condition" for the delta-function, 



lim 

e-S-0 



dz 



2 „ 2 c 2 ±c-6 4 . 
-p 2 - «9 2 2 + 5 + -d z 



F ± (p, z, z') = i(kz 



/\4 



(115) 



(116) 



provide two relations. Two additional conditions that are required for an unambiguous 
determination of the propagator come from the boundary conditions on the branes. One 
has to enforce that the fermion propagator only contains the zero-mode of the correct 
handedness. These conditions are most conveniently found by explicitly writing (IllOp in 
matrix form. For the propagator of the SU(2) singlet field E one finds 



A(p,z,z') 



{d z -\- f)F-(p, z, z')t -^p,F+(p, z, z') 
-(7%i^(p, z, z') (-d z + f - l)F+(p, z, z')t 



;il7) 



The off-diagonal elements of (I117P conserve the handedness during propagation; the 
diagonal elements lead to mixing of modes of different handedness. Since in the mode 
representation (see (I123aj) below) the upper-right submatrix corresponds to the pure 
"left" case, its entries must vanish on the branes for all values of as follows from ( I88p : 
thus 



F E (P> Z ') 



z=l/k,l/T 



0. 



(118) 



Similarly, the upper-left submatrix must vanish for z = l/T,l/k as the left-handed 
components of the singlet field must vanish!^! This gives the second boundary condition: 



[d z -- z -- z }F-(p,z,z') 



0. 



(119) 



z=l/k,X/T 

Analogously, the boundary conditions for the doublet fermion propagators are given by 



F L (p,z,z') 



[-d. + ---]F+(p,z,2>) 

z z 



z=l/fc,l/T 



z=X/k,l/T 



o, 

0. 



(120a) 
(120b) 
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Note that this is not true for z' = X/T, 1/k. 



45 



It is now straightforward but tedious to determine the coefficient functions A and B 
which again can be expressed in terms of Bessel functions J and Y . We give the results 
in a more concise form below. 

For our calculations it is most convenient to decompose the fermion propagators into 
their RL, LR, LL and RR chiral components: 

x,y) =- F~(p, x, y)fP L - F£(j>, x, y)fP R 

+ d + Fx(p, x, y)P L + d'F+ip, x, y)P R , (121) 

where the X can be either E or L. The differential operators 

d+ = d x - - - - d ~ = -d x + ---. (122) 

always act on the first bulk position argument (x in (j!2ip ). The decomposition (11211) 
also allows for a simple representation in terms of Kaluza-Klein mode sums: 

F£(p,x,y) 

d~F+(p,x,y) 

d + F^(p,x,y) 

Analogous expressions hold for the singlet field E. 

The chiral components d ± F T and F^ 1 are scalar functions; the Dirac and Lorentz 
structures of the propagator are contained only in the projectors and f. Since the final 
evaluation of the loop integrals is performed numerically we need the analytic expressions 
for these scalar functions after a Wick rotation has been performed. In our matching 
calculation this rotation must be performed after the (Minkowski) expressions for the 
Wilson coefficients have been extracted and the Dirac and Lorentz algebra are already 
done. The Wilson coefficients are then given as integrals over scalar functions. The Wick 
rotation effectively amounts to the replacement p 2 — > —p 2 or p — > ip with p = \/p 2 , where 
p is the loop momentum. The Bessel functions J and Y are replaced by the modified 
Bessel functions I and K of the first and second kind, respectively, via 

J u (ix) = i v Ux) Y u (ix) = i v+l I u (x) - - r v K v {x) . (124) 

To write the chiral components of the fermion propagator in a compact way we make 



= T,ff ) ^^^ft ) (y) (123a) 

n f n 

= H^^)-f^^{y) (123b) 

n n 

= J2gP(x)^^fP(y) (i23c) 

n 

-Ef'^J^^Hy) (123d) 
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use of the following functions which were first introduced in [15]: 

S+(p, x, y, c) = I c+1/2 (px)K c+1/2 (py) - K c+1/2 (px)I c+1/2 (py) (125a) 

S-(p, x, y, c) = I c -i/2(px)K c _ 1/2 (py) - K c _ 1/2 (px)I c _ 1/2 (py) (125b) 

S+(p, x, y, c) = I c+1/2 (px)K c „ 1/2 (py) + K c+1/2 (px)I c -i /2 (py) (125c) 

S-(p,x,y,c) = I c -i/2(px)K c+1/2 (py) + K c _ 1/2 (px)I c+1/2 (py) . (125d) 

Additionally, the following functions are useful when Taylor expansions of the propaga- 
tors are needed 



S'+(p,x,y,c) =^^S + (p,x,y,c) - ^S^(p,x,y,c) + ^S+(p,x,y,c) (126a) 

S'_(p,x,y,c) = 1 ^ 2 2C S-(p,x,y,c) - ^S + (p,x,y,c) + ^S-(p,x,y,c) (126b) 

S+(p,x,y,c) =^S + (p,x,y,c) - ^S^(p,x,y,c) + ^S + (p,x,y,c) (126c) 

S'_(p,x,y,c) =^S_(p,x,y,c) - ^S + (p,x,y,c) + ^S^(p,x,y,c) . (126d) 



The corresponding expressions in Minkowski space can be recovered with help of the 
following substitution rules: 

S + [I,K] ~S+[J,Y] ~S + [Iu,K,} -> ip + [J„, {-IY- C+1 ' 2 Y,} (127a) 

S_[I,K] ~S-[J,Y] S-.[I V ,KA -> -i^S^i-iy-*- 1 ^] (127b) 

where the brackets [•, •] show how I and K have to be replaced with the Bessel functions 
J and Y, respectively. 

For the right-handed singlet fermion field the scalar functions are given by 

„+, s rv tk*x 5 / 2 y 5 / 2 S + (p, x, l/T, c E )S + (p, y 1/fc, c E ) 
(P ' *' y) = " e( * " V) S + (p,l/T,l/k,c E ) 

Q , \ ^V /2 y 5/2 g + (p, y, l/T, c g )£ + (p, x, 1/fc, eg) , 

— &{y — x — — , , r ^izaaj 

S+(p, l/T, 1/fc, cs) 

, zfcV/ 2 y 5 / 2 g_(p, x, l/T, c g )g_(p, y, 1/fc, eg) 
(P ' X > y) = Q(X ~ V) S + ( P A/T,l/k,c E ) 

S+(p,l/T,l/k,CE) 

itfxWyWS+fax, l/T,c E )S^p,y, l/k,c E ) 
d + F E (p, x, y) =p0s- y) TaFTTI: ^ 
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^ n( ik*g*/*yS/*S- (p, y, 1/T, c E )S + (p, x 1/k, gg) 
+ pB(?/-x) . 1/^1/1 \ ( 12 8c) 



d F^(p,x,y) = pQ(x - y) 



S-(p,l/T,l/k,c E ) 

ik^/ 2 yV 2 ~S_(p,x, l/T,c E )S + (p,y, l/k,c E ) 
S+(p,l/T,l/k,c E ) 



^ n( ik*xV 2 yWS + (p, y, 1/T, c E )S.(p, x, 1/fc, eg) 

+ p6(y - x) — r . 128d 

S+(p,l/T,l/k,c E ) 

For the left-handed doublet field L we find 

itfxWyWS+fa x, 1/T, c L )S + (p, y, 1/k, c L ) 



Ft(p,x,y) = Q(x-y) 

+ e(y-x) 

F£(p,x,y) =-B(x-y) 
- Q(y - x) 



S^p,l/T,l/k,c L ) 

ik A x^l 2 ~S + {p, y, 1/T, c L )S + (p, x, 1/k, c L ) 
S4p,l/T,l/k,c L ) 

ik A x^ 2 y^ 2 S^p,x, l/T,c L )S.(p,y, l/k,c L ) 
S-(p,l/T,l/k,c L ) 

itfxBfiyWS-fa y, 1/T, c L )S.(p, x, 1/k, c L ) 



d + F L (p, x ,y) = - pQ(x - y) 

— pQ(y — x) 
d ~F^(p, x, y) = - pQ(x - y) 

— pQ(y — x) 



S.(p,l/T,l/k,c L ) 
ik*x*/*y*/*S + (p, x, 1/T, c x )£L(p, y, 1/k, c L ) 



5_(p, 1/T, 1/k, c L ) 

itfxWyWS-jp, y, 1/T, c L )S + (p, x, 1/k, c L ) 
S4p,l/T,l/k,c L ) 

lk ^/2 y 5/2 S _ ^ x? 1/T> CL )s + (p, y, 1/k, c L ) 

S-(p,l/T,l/k,c L ) 

ik^' 2 y^ 2 ~S + {p, y, 1/T, c L )S.(p, x, 1/k, c L ) 
S4p,l/T,l/k,c L ) 



;i29a) 



(129b) 



129c) 



(129d) 



Note that here p = y p 2 = yp\ + p\ + p\ + p\ is the Euclidean momentum. 
A. 3. 2 Gauge boson propagators 

The derivation of the propagators for the gauge sector follows [9]. The gauge is used 
to disentangle the orbifolding-even vector components of the 5D gauge field from the 
orbifolding-odd fifth component. The propagator for the "vector part" is then determined 
by the equation 



A. pu (p, z, z) = ikz 5(z — z) 



(130) 
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Substituting the ansatz 



A^(p,x,y,£) = A ± (p,x,y) (r) 



P 



+ — 5-A||(p,x,2/,0 



(131) 



into ( 11301) gives two equations 

(p 2 + z9^-9 z )A_l(p, z, z) = —ikz 8(z - z) 



(132a) 



(-p 2 - zd z ^d z )A ± (p,z,z') - 1 Up 2 -£z<9,~£U A,|(p,z,/,0 = 0. (132b) 

Putting A||(p, z, z', £) = A±(p/y/£,z,z') the second equation is automatically satisfied, 
if A_l(p, z, z') solves the first one. Eq. f |132aj) can be solved with the same strategy 
employed before for the fermion case. Here the boundary conditions read 



d z A ± (p,z,z')\ z=1/kA/T = 0. 



(133) 



The propagator of the fifth component of the gauge field is denoted by A 5 (p, x, y, £) 
and follows readily from 



p 2 + £d z zd z 



A 5 (p, z, z) = ikz 8(z — z) 



and the boundary condition 



A 5 (p,z } z')\ z=1/kt 



l/T 



0. 



(134) 



(135) 



Again we only give the Euclidean propagators: 

ikxyS + (p,x, l/T, l/2)S+(p,y, 1/k, 1/2) 
Ax(p, x, y) = e(x - y) s^PA/T, 1/k, 1/2) 

_ v *A;xy,g + (p, y, l/T, l/2)S + (p, x, 1/k, 1/2) 
+ [V X) S-(p } l/T, 1/k, 1/2) 

Ay (p, x, y, = A ± (p/ v^, x, y) 

*A;:ryS_ (p, x, l/T, l/2)5_(p, y, 1 A, 1/2) 



(136) 
(137) 



£A 5 (p,x,2/,f) = 9(x - y) 



+ 0(y - x\ 



S-(p, l/T, 1/k, 1/2) 

ikxyS.jp, y, l/T, l/2)S_(p, a, 1/fc, 1/2) 
S-(p,l/T,l/k,l/2) 



;i38) 



where p = p/ The Minkowski space expressions can be recovered via (I127p . The SM 
contribution to g — 2 is, from point of view of the matching calculation at the KK scale, 
a long-distance contribution. As described in the main text, this contribution is most 
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effectively removed by subtracting the contribution of the SM gauge boson from each 
5D boson propagator, which is equivalent to subtracting the zero mode. We therefore 
define the zero-mode subtracted (ZMS) propagators 

Ai MS (p,x,y) = A ± (p,x,y) -±fW( x )ff)( y ) (139) 

Af MS (p,x,y,£) = A™ s (p/^,x,y) . (140) 

The boundary conditions for the fifth component of the gauge fields prevent the existence 
of a massless mode; no subtraction is necessary for A 5 (p, x, y, £). 

A. 3. 3 Higgs propagator 

The Higgs field is IR brane localized. In the unbroken phase the Higgs propagator has 
a tachyonic mass /j, 2 . Since fi 2 ~ v 2 <C T 2 we can treat the mass term as a perturbation 
of order v 2 in the calculation of the matching coefficients, and work with a standard 
massless scalar 4D propagator for the Higgs, 

A„(z, „)„= J |^ e -*->^= (141) 

where m, n = 1,2 are SU(2) indices. 



B Expanded propagators 

If the virtuality of the loop momentum is of the order of the KK scale T one can expand 
the propagators in powers of the small external momenta. We effectively perform a 
Taylor expansion of the form 

d 

A(p - k, x, y) « A(k, x,y)-2p-k ^A(^, x,y) + ... (142) 

where all momenta are still Minkowskian and A can be any propagator. Using the identi- 
ties for derivatives of Bessel functions we obtain analytic expressions for the derivatives. 
In the following we denote -M^A{k,x,y) by A{k,x,y). All expressions are only valid 
after the Wick rotation has been performed. The expressions for fermion propagators 
read@ 



ik A x 5 / 2 y^ 2 Q{x - y) 
S + (p,l/T,l/k,c) 

15 To shorten the notation, we denote the bulk mass parameter by c, where c = ce in the SU(2)-singlet 
fermion propagators and c — cl for the doublet ones, as in (|128p . (|129p . 



S' + (p, x, 1/T, c)S+(p, y, l/k, c) + S+(p, x, 1/T, c)S' + (p, y, 1/k, c) 
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S+jp, x, 1/T, c)S+(p, y, 1/fc, c)S' + (p, 1/T, 1/fc, c) 
S+(p, 1/T, 1/fc, c) 



S'+(p, y, c)S+(p, x, 1/fc, c) + S+(p, y, 1/T, c)S' + (p, x, 1/fc, c) 

(143a) 



1/T, 1/fc, c) 

S+(p, y, 1/T, c)S+(p, x, l/k, c)S' + (p, 1/T, 1/fc, c) 



S+(p,l/T,l/k,c) 



F E (p, x, y) = 
ik*x 5 / 2 y 5/2 ®(x - y) 



S'_(p,x, 1/T, c)S-(p,y,l/k,c) + S_(p,x, 1/T, c)S'_(p,y,l/k,c) 



S + (p, 1/T, 1/fc, c) 
S-(p,x, 1/T, c)S-(p,y, l/k,c)S'(p, 1/T, 1/fc, c) 



S+(p, 1/T, 1/fc, c) 

+ ^^17^/^ (^-(P. I/» VT c)5_(p, 1/fc, c) + S_(p, y, 1/T, c)^(p, x, 1/fc, c) 
S_(p, y, 1/T, c)S-(p, x, 1/fc, c)5;(p, 1/T, 1/fc, c) 



S + (p, 1/T, 1/fc, c) 



(143b) 



T L (p, x, y) = 

^4^5/2^5/20^ _ y j 



S'_(p,x, 1/T, c)S_(p,y, l/fc,c) + S_(p,x, 1/T, c)S'_(p,y, 1/fc, c) 



S_(p, 1/T, 1/fc, c) 

S_(p,x, l/T,c)S-(j>,y, l/k,c)S'_(j>, 1/T, l/fc,c) 
S_(p, 1/T, 1/fc, c) 

ifcV/V/ 2 0(y - x 



(p, y, 1/T, c)S_(p, x, 1/fc, c) + S_(p, y, 1/T, c)#_(p, x, 1/ fc, c) 

(143c) 



S_(p, 1/T, 1/fc, c) 
S_ (p, y, 1/T, c)S_(p, x, 1/fc, c)tf_(p, 1/T, 1/fc, c) 



S_(p, 1/T, 1/fc, c) 



F£(p,x,y) = 
ik*x 5 / 2 y 5 / 2 e(x - y) 



S'+(px, 1/T, c)S+(p, y, 1/fc, c) + S + (px, 1/T, c)5+(p, y, 1/fc, c) 



S_(p, 1/T, 1/fc, c) 

_ S + (px, 1/T, c)S+(p, y, 1/fc, c)S'_(p, 1/T, 1/fc, c) 
S_(p, 1/T, 1/fc, c) 

+ %^^Ti/kc) (^ +(P ' y ' VT ' C) ^ +(P ' 1/A; ' C) + y ' 1/T ' C) ^ (P ' C) 
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S+jp, y, l/T, c)S+(p, x, 1/fc, c)S'_(p, 1/T, 1/fc, c) . , . 

5_(p,i/r,i/fc, c ) 1 ( -- )<lj 



d~F+(p,x,y) = 

' lV ' V ! B( ''' " _ ^ ^-(p,*, l/V.H.S. (/,//• 1 //■••'•) + 5_(p,x, l V.rl.V. (/,//. l/fc,c) 



S+(p,l/T,l/k,c) 

S_(p, x, 1/T, c)S + (p, y, 1/fc, c)S;(p, 1/T, 1/fc, cj 
5 + (p, 1/T, 1/fc, c) 

+ %P ^f 2 f!y}n~^ (S' + (P, V, 1/T, c)5_(p, 1/fc, c) + S + (p, y, 1/T, c )#_(p, x, 1/fc, c) 



5 + (p, 1/T, 1/fc, c 

s+(p, V, 1 / T , c)S-(p, x, 1 A, c)s+(p, 1/T, 1 A? c ) 
S+(p,l/T,l/k,c) 

1 

Yp 



2 d F+(p,x,y) (144a) 



d+F E (p, y) = 

""■' / ' ' ' !j) f S'Jp, x, 1/T, c)S_(p, y, 1 IK c) + S + (p, x, l/T,c)S'_( P ,y,l/k,c) 



S+(p, 1/T, 1/fc, c) 

S + (p, x, 1/T, c)S4p, y, 1/fc, c)S' + (p, 1/T, 1/fc, c) 
S+(p,l/T,l/k,c) 



5 + (p, 1/T, 1/fc, c) 

S-(p, y, 1/, c)g+(p, x, 1/fc, c)S' + (p, 1/T, 1/fc, c) 
S + (p, 1/T, 1/fc, c) 

1 

27 2 



-— d+i^(p,x,y) (144b) 



d + F L (p, x, y) = 

' y) ^S' + (p,x, l/T,c)S-(p,y, l/kc) + S+(p,x, l/T,c)S'_(p,y, 1/fcc) 



5_(p, 1/T, 1/fc, c) 

1/T, c)g_(p, y, 1/fc, c)fl_(p, 1/T, 1/fc, c) 
5_(p, 1/T, 1/fc, c) 
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>l ' k ' ' " ' " H{ !l " ' ' 1 { S'_{p, y, l/T, c)S + (p, x, l/k, c) + S. (p, y, 1/T, c)S> + (p, x, 1/k, c) 



S-(p,l/T,l/k,c) 

S_(p,y, l/T, c)S + (p,x, l/k,c)S'_{p, l/T, l/k,c) 
S-{p,l/T,l/k,c) 

1 

2p 



—d + F^(p,x,y) (144c) 



d F£(p,x,y) = 

^ ' ' ' ~ U) f SL( P ,x, l/T,c)S + (p,y, l/k,c) + S-(p,x,l/T,c)S' + (p,y,l/k,c) 



S-(p,l/T,l/k,c) 

S-(p, x, l/T, c)S + (p, y, l/k, c)S'_(p, l/T, 1/k, c 
S-(p,l/T,l/k,c) 

ik 4 x 5/2 y 5/2 Q(y - x 



S+{p, y, l/T, c)S-(p, x, l/k, c) + S+(p, y, l/T, c)]S'_{p, x, l/k, c) 



S-(p,l/T,l/k,c) 

S + (p, y, l/T, c)S-(p, x, l/k, c)SL (p, l/T, l/k, c) 
S4p,l/T,l/k,c) 

1 

2p 



— d-F+{p,x,y) (144d) 



For the gauge boson propagators we need: 
Aj_(p, x, y) = ikxy Q(x — y) 

' S' + (p, x, l/T, l/2)S + (l, y, l/k, 1/2) S + (p, x, l/T, l/2)S' + (l, y, l/k, 1/2) 
SL(p, l/T, l/k, 1/2) S-(p, l/T, l/k, 1/2) 

S + (p, x, l/T, l/2)S + (l, y, l/k, l/2)S'_(p, l/T, l/k, 1/2)' 



S_(p, l/T, l/k, 1/2)2 
+ {x^y} (145) 



A 5 (p,x,y)= t -^-Q(x-y) 



' S'_(p,x, l/T, 1/2) S_(p,y, l/k, 1/2) Sjj^x, l/T, l/2)S'_{p,y, l/k, 1/2) 
SL(p, l/T, l/k, 1/2) + S4p, l/T, l/k, 1/2) 

S-(p, x, l/T, l/2)S_(p, y, 1/fc, l/2)ff_(p, l/T, l/k, 1/2) \ 
5_(p, 1/T, 1/2)2 J 



53 



+ {x y} 



(146) 



C Explicit expressions for the diagrams 

In this section we provide explicit expressions for the diagrams relevant to the matching 
of the dipole operator. To simplify the expressions, we consider the diagrams with an 
external photon, that is the linear combination cwB^ + sw\W^'^ of diagrams with an 
external hypercharge and SU(2) gauge boson. We furthermore assume that the external 
Higgs field (represented by a grey square) has a vanishing upper isospin component, that 
is we pick up the lower component only, since this is relevant when <3> is replaced by its 
vacuum expectation value. These simplifications allow us to drop diagrams, which are 
proportional to cwg$ — s w gs = 0. Using e$ = cwg'$ = s w g^ it is always possible to extract 
the 5D electromagnetic coupling. The electric charge of a fermion ip in units of positron 
charge is denoted by Q^. In our convention the hypercharges are given by Yl — — 1, 
Y e = —2, Y$ = 1, and = + W22 = — 1- We use the short-hand p = p — I and 
ft = p' — I for momenta in the loops. 

The labels of the diagrams correspond to those given in figure [21 The first expression 
for each diagram does not make use of any simplification. Vertex factors and propagators 
are written in the form as given in appendix |X] and no terms have been omitted. The 
second expression already includes some simplifications. In particular, the terms that 
vanish due to the presence of chiral projectors after multiplying out the chiral components 
of the fermion propagators have been dropped. This second representation is exact only 
if the diagram does not contain an insertion of the Higgs field into an external line, 
since we applied the replacement (1471) . The (small, see discussion in the main text) "off- 
shell" terms present in these diagrams is not shown explicitly and must be added. These 
contributions can be readily obtained from the unsimplified expressions by replacing the 
external fermion propagator by the zero-mode propagator only, e.g., in B2a one would 
replace 

Af(p,x,l/T) f^(x)p^(l/T). (147) 

One then applies standard Dirac algebra and simplifications as for the other term shown 
explicitly, including an expansion in the external momenta p, p' to the appropriate order. 
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c 'W n ^ fz!^ (5D) XlXe. f 1/T dx f 1/T d JL f 1/T dz f dH 

iWs) y k s ) Vii 2 2 J yk {kx y J i/k {ky y ^ {kz y J (27r) 4 

f^(z)ff(y)g^(x)e^AZ s (l,x,z) 

xLi(j/)P R7p Af(p' - Z, z, t/) 7tl Af(p - /, y, l/T)Af (p - Z, 1/T, x^PrE^p) 
g'iesQ^YEy^T 3 dx f^ T dy f 1 ^ dz f dH 



4fc3 ; 1/Jfe (A;x) 4 (%) 4 Vi/fc (kz) 4 J (2tt) 4 

xLiCp 7 ) [Flip 1 , z, y)F+{p 1 y, l/T)F^(p, 1/T, x) {i P (p"-?h^} (p - Z) 2 + 
rf + F L " (p', 2, y)d-F+(p, y, 1/T)Fe.(p, 1/T, a;) { 7p7/1 (^-/) 7l/ }] P^p) (148) 

r 'W n \ (~ lT ^ ( 5 °) dx f 1/T dy f 1/T dz !J-L 

W y k3 j y l3 2 2 (%)4 ( ^ )4 y (27r)4 

(y)^ ) (x)e^A£i B (Z,x l z) 
xLMP^Afip' - I, z, l/T)Af(p' - /, 1/T, y) lfM Af (p - Z, y, x^PrE^p) 

&e s Q li Y L Y 1 ty$ D) 7* dx dy f 1 ^ dz f dH 



4F (fc^ A/* (%) 4 A/* (^) 4 7 (2vr) 4 

xL,(p') [f+(p',z, l/T)F^.(p', l/T,y)F^{p,y,x) {iMp'-fM (p' " 0" + 
T+(p', z, l/T)d~F+ (p', 1/T, y)d+Fg(p, y, x) { lp (f-/)^}] P R E 3 (p) (149) 
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R9 ,. , w . n , f~iT 3 \ m Y L Y L [ l ' T dx f^dy [ 1/T dz f dH 



xLMPn^Afo' - I, z, y) 7 ,A i L (p - Z, y, x) 7l/ Af (p, x, l/T)P R Ej(p) 



r 1/T dx 


f^dy 


f 1 ^ dz 


f dH 


L/k {kxY J 


i/k ( W J 


i/ k (kzy J 


(27T) 4 



Ak 3 

fff(z)f?\y)g^(l/Ty»AZs(l,x,z) 

x L t (p') [F+ (p\ z, y)d + FZ i (p, y, (p, x, 1/T) { Tp (/-/) 7^} + 

rf + F L ~ (p', 2 , y)F£ i (p, y, (p, a;, 1/T) {i P i^-^} ] P^- (p) (150) 

R9 , / W - n x f~ iT3 \ (sd)Y e Y e f l ' T dx [ l ' T dy f l ' T dz f dH 

B2b = (ig 5 ) (,e 5 g,) ^— j y < TT / /fc / 1/fe (^ / (2^ 

xLi(jf)P R Af(p', 1/T, z) lt Af{p' - I, z, y)^Af(p - I, y, x^PrE^p) 
&e*Q*YhT )l * f 1/T dx f l/T d y f 1/T dz f d ' 1 



J 1/k (kxr Ji/k (kyf J 1/k (kz)* J (2vr)4 
f^(l/T)f\ \y)gP j (x)e^AZ s (l,x } z) 

xLi{p) d~F+ (p\ 1/T, z)d + F^ (p, z, y)F E ~. (p, y, x) {l P l^-/hu} + 

d~F+ (p\ 1/T, z)F+ (p', 2, y)d + F^ (p, y, x) {7 p (p"-/)7m7,} 1 PrE 3 (p) (151) 
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iT 3 \ ( 5 D) 



Yl Ye f 1/T dx f 1/T dy f 1/T dz 



dH 



{9 ' 5)2{ie5Q ^ ' *3 )*« 2 2 J 1/k {kxf J 1/k (%) 4 J 1/k {kzY J (2tt)< 
f^(z)f^(y)g^(x)e^A 5 (l,x,z) 

xL t (p')P Rl5 Af(p' - I, z, yh,Af(p - I, y, l/T)Af (p - I, 1/T, x)y 5 P R Ej(p) 



4k 3 



r 1/T dx 


f^dy 


f 1 ^ dz 


f dH 


h/k (k x y J 


i/k (%) 4 J 


i/k ^zY J 


(2vr) 4 



f^(z)f^(y)g^(x)e^A 5 (l,x,z) 

xLi(f) [d-F+tf,z,y)F%fay, l/T)cTT+ (p, 1/T,x) {7^-/)} + 
F7(p',z,y)d-F+(p,y, lffltTF+M 1/T,x) {(/-/h^l 1 PsPjip) 



(152) 



Y L Y E f 1/T dx f 1/T dy f 1/T dz 



dH 



k 3 ; 2 2 J 1/k (kxf J 1/k {kyY h/ k (kzY J (2tt)* 

fS(z)f!?\y)g ( £(x)e*»A 5 (l,x,z) 



xL^Pja^ip' - l,z, l/T)Af(p' - I, l/r,y) 7 „Af (P " l,V,x) l5 P R Ej{p) 



Ak 3 



r 1/T dx 


f^dy 


f^ T dz 


f d 4 l 


h/k (kx) 4 J 


i/k (kyY J 


i/fc (kzY J 





f^(z)ff(y)gP j (x)e^A 5 (l,x,z) 

xLtip') [d~F+ {p\ z, l/T)d-F+ (p', 1/T, y )F+ (p, y, x) {7 M (*-/0} + 

d~F+ (p', z, 1/T)F£ (p', 1/T, y)rf-T+ (p, y, x) {{p"-/h»} } P R E s {j>) (153) 
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R/1 < n , iSD )Y L Y L f l ' T dx f 1 ^ dy [ l ' T dz f dH 

/ffto/f (^(l/T^^,^) 

xLMPr^A^p' - I, z, yh,Af(p - I, y, x) 75 Af(p, x, l/T)P R Ej{p) 



f 1 '* dx 








L/k {kxY J 






(27T) 4 



4A; 3 

/1?(*)4 0) (^l/T)e*<*A 5 (Z,^) 

xL.(p') [eTT+(p', 2, 2/)F+(p, y, x)d-F+(p, x, 1/T) {7^-/)} + 

F L -(p', z, y)d-F+(p, y, x)d-F+(p, x, 1/T) {{p"-f) lfl } ] P B ^(p) (154) 

-tT 3 \ (5D) Y E Y E f 1/T dx f 1/T dy f 1/T dz f dH 



B4b = (, 5 ) K g M ) ^— j 4 YY L L (35? ./ 



2 2 y i/fc (fca;) 4 y i/fc (%) 4 7 1/fc (A^) 4 7 (2vr) 4 
f^(l/T)ff\y)g^(x)e^A 5 (l,x,z) 

xLiWPnAfip', 1/T, *) 7s Af (p' - I, z, y)^Af(p - I, y, x^PrE^p) 
g'iesQ.Ylyf^ dx f x ' T dy f^ T dz f dH 



4fc» J 1/k (kx)^ Ji /k (%) 4 Ji/k (kzY J (2tt)4 

f^{llT)f^\y)g^{x)e^A^S^z) 

xLiip') [d-F+ip 1 , l/T,z)d-F+(p',z,y)F+(p,y,x) { 7 „ (*-/)} + 

<TF+ (p', 1/T, z)Fg (p', z, y)d-F+ (p, y, x) {(/-/) 7 J } PrE^p) (155) 
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B5a = (ig' 5 Y(ie 5 Q 



2 2 



1T 3 



(5D) 

Vij 



X ' T dy f 1 ^ dz 



dH 



k 3 J ^ J 1/k {kyY Ji/k (kz)* J (2nY 
/£ ) («)4 o) (y)^ ) (l/T)^A H (0A22 B (i I l/r,2 ! ) 
xLi(jpf)Pja a Af(p' - I, z, y)^Af(p - I, y, l/T)(-/ /3 )P^(p) 
g'ie.Q.Y^yf^T 3 [V T dy [V T dz f dH 



4^ 3 h /h {kyY J 1/k (kzY J (2tt) 4 

/£ ) W4 0) (2/)4 ) (l/T)e^A H (/)A^ s (U/r^) 
xLi(p') [d+F^', z, y)d~F+(p, y, 1/T) ( 7q7m ) l p + 
Fl{p', z, y)F+(p, y, 1/T) {^(p"-/)^-/)} / J P R ^.(p) 



B5b = (ig' 5 Y(ie 5 Q 



YeY* 
2 2 



iT 3 



A- 3 



(5D) 



l ' T dx 



1/T 



dy 



f^^l/T)f\°\y)g^(x)e^A H (l)A^l J x 1 l/T) 



xLiip^PnAfip' - 1, 1/T, y) ltl Af (p - Z, y, x^PrE^p) 



g^Q^Yzyf^T 3 dx 



1/T 



d 4 / 



4fc3 y i/fc (A;x) 4 A/ fe (%) 4 7 (2vr) 4 

/ff(VT)4 0) (y)^(x)e*"A ff (OA£ B (i,x,l/r) 
xZ^p') \d~F+ (p', 1/T, y)*^. (p, y, x) ( 7/i7a ) ^ + 
P£(p% l/T,y)F^(p,y,x) {(/-/) 7/ ,(^-/) 7a } J P^p) 



(156) 



i/A (M 4 A/* (%) 4 7 (2vr) 



d 4 / 



(157) 



59 



The amplitudes for the two diagrams above, Wl and W2, can be obtained from the 
amplitudes for diagrams B2a and B4a, respectively, by replacing g' 5 2 Yj* with g\ in ( |150p 
and f TTMjl . 




The amplitude for the diagram above, W3, can be obtained from the amplitude for 
diagram B5a by replacing g'^Y^Y L with g\ in (I156p . 
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W4 = (^ 5 ) 2 (-e 5 ) 



-iT c 



(5D) 

Vij 



22 



l ' T dx 



1/T 



dy f 1/T dz f d 4 l 



i/k (kx) 4 J 1/k (ky) J 1/k (kz) 4 J (2tt) 



{(p' - p) - (p - 0)vty*A + {(P - - (Z - p')}fiVu\ 

+ {{l-p') - (P'-P)hvj * L,(p / )Pi?7«Af(/,z,x) 7/3 Af (p,x, l/T)P R Ej(p) 



95 e 5V 



IJ 



r 1 ^ dx 


ray 


( l ' T dz 


f dH 


li/k (kxy J 


i/fc (ky) J 


i/fc (kz) 4 J 


(2vr) 4 



2k 3 

/i?W4 0) (y)g^(l/T)e^A^ s (p,x,y)AZ s tf,y,z) 

(p' - 2p + Z)^A + (p - 2/ + p%^A + (/ - 2p' + p)a^ 
xA(p') k^" (/, z, v)d-F+(p, x, 1/T) { 7a7/3 } 1 P R E,(p) 



W5 = Gfc) 2 (e 5 ) 



iT 3 



k 3 



(5D) 



2 2 



22 



(kx) 4 



i/fc 



1/T _^ 
i/fc (%) 



l ' T dz 
i/fc 



/i° ) (^)4 0) (|/)4 ) (l/r)e^A 5 (p - l,x,y)A 5 ( P ' - l,y, z) 
x{(p-l)-(l-p')}„ 

xL i (p / )P R 7 5 Af(/,2:,x)7 5 Af(p,x,l/T)P R ^(p) 



2k 3 



[V T dx 




f 1 ^ dz 


f d 4 l 


A/fc (kx) 4 J 


i/fc (ky) J 


i/fc (kz) 4 J 


(2tt) 4 



f { S(z)f^(y)gf.(l/T)e^A,(p - l,x,y)A 5 (p' - l,y,z) 

x (p - 21 + p% L 4 (p') [<TF+(Z, z, x)d-F+(p, x, 1/T)] P R Ej(p) 



(158) 



d 4 l 



(159) 
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W6a = (ig 5 )g5(-ie 5 ) 

KO)/ s (o) 



-iT 3 



k 3 I 3 



,(5D) 



22 



^dz /-^dy Z^cfe 



V* (kxY J 1/k (ky) J 1/k {kzY J (2tt) 



/i7W^(l/T)e^A 5 (p-U,y) 



A&fo' - l,y,z)^ff(y) - f^\y)^KL{p' - l,V,z) 



xLiip^PsrfaAf (I, z, x) 75 Af (p, x, llT)P R Ej{p) 



2k 3 



[ l ' T dx 


f^dy 


fV T dz r 


ll/k (kxY j 


i/fe (ky) J 


i/k (kzY J 



d 4 l 



fS(z)fi°\y)gP j (l/T)e;A 5 (p-l,x,y)x 
xU(p') \f+(1, z, x)d~F+(p, x, 1/T) ( 7a /) P R Ej(p) 



^-AZ s (p'-l,y,z) 



W6b = gS9b){-i^) 



iT 3 



k 3 



(5D) 

Vij 



r a r b 



^dx /^dy /-^^ 

l/k 



22 A/fc (^) 4 A/fc (ky) Jx/h (kz) 4 J (2tt 



f^(z)g { S(l/T)e^A 5 (p'-l,y,z) 
XV fa 

xLi^PflTgAf (/. - . .»•)-„, Af(,>. .r. l/T)P H Ej(p) 



A^ s (p - l,x,y)^\y) - f^(y)^-A^ s (p-l,x,y) 



dH 



2k 3 



r 1/T dx 


f^dy 


f 1 ^ dz f 


li/k (kxY J 


i/k (ky) J 


i/k (kz'Y J 



f^(z)ff\y)g^(l/T)e;A 5 (p'-l,y,z)x 
xLiip') \F Li (l, z, x)d-F+(p, x, 1/T) 7a ) P R Ej(p) 



^-AZ s (p-l,x,y) 



dH 



(160) 



dH 



(161) 
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W7 = ('ijfc) 2 (e 5 ) t 



Y -± T 1 + V 
2 2 4 



22 



iT 3 



A- 3 



(5D) 

Vij 



l ' T dz 



l/k {kzY J (27T 



f^z)f^(l/T)g^(l/r)^A H (p - QA£b(p' - J, 1/T, z) 
xL i (j/)P R 'y v At{l 1 z 1 l/T)P R E j (j>) 
dz f a 11 



2k 3 J 1/k {kzY J (2tt) 4 

/i! ) (^)4 0) (i/T)4 ) (i/T) e ;A / ,( P - 1)AZ S (P' - i, *) 

x£,(?/) z, 1/T) ( 7 ,/) 1 P^(p) 



d 4 / 



(162) 



W8 = (ig 5 ) 2 (ie 5 ) - 



2 + 2 



~2 



iT 3 



A- 3 



(5D) 

Vij 



l ' T dz 



i/k 



dH 



{kzf J (2tt) 4 

f { S{z)ff\l/T)g%{l/T)e^A H {p - l)A H (p> - l)A^{p> - 1, 1/T, 
xLifflPma&til, z, 1/T) (I - p')p {(l-p) + (l- p')} u P R E 3 (p) 



2A; 3 



1/k {kzY J (2vr) 



dH 



f^{z)ff\l/T)g%{l/T)e^A H {p - l)A H {p' - l)A^ s (p' ~ h *) 
xL t {p') \F+{l,z, 1/T) ( 7q /) (p' - l) p (p + p>- 2l)]p R E j (p) (163) 
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yy 



ab3 



-iT 3 



(5D) 

Vij 



1/T 



dy f 1/T dz f d 4 l 



K) 2 ( e 5) ^ 2 2 j ^ . y k3 j yiJ {ky) {kz) , j (2?r)4 

f { £(z)f?\y)g®(l/T)e*»A H (p - l)A^ 3 (p - 1, 1/T, y)A™ s (p' - I, V, z) 

{(p' -p)-(p- 1)}vVix\ + {(P~ l )-( l - p'))hVpX 
+ {{l-p')- ip' - p)}xvJ x W)P Rla Af(l, z, 1/T)(l - p)pP R E j (p) 



X 



9ie 5 yf ] T 3 f^dy dz 



d 4 l 



2k 3 J 1/k (ky) J 1/k {kzY J (2tt) 4 
f^(z)ff\y)g^(l/T)e^A H (p - l)A^ s (p ~ I, l/T,y)AZ s (p' ~ l,y,z 

(p' -2p + l)^^ + {p-2l + p')^i] uX + (l- 2p' + p) x ri IMU 
xL t (p') \F+(l, z, 1/T) ( 7o /) (I - p) p ] PrE^p) 



x 



(164) 



(fl , 5)(^5)(-ie 5 ) 



iT 3 



k 3 



(5D) 

Vij 



yy 



J 22 



(ky) J 1/k {kzY J (2tt) 



ci 4 / 



xL^Pr^A^z, l/T)(l-p) p P R E j (p) 

gjenyjfT 3 fV T dy_ f l ' T dz f dH 
2k 3 J 1/h (ky) J l/k (kz) 4 J (2tt) 4 

f i Lhz)ff\y)gf j (l/T)^A H (p - l)A h (p' -l,y,z 




x 



dy A^ s (p-l,l/T,y) 



xL t (p') d-F+(l,z,l/T)(l-p)p P R E 3 (p) (165) 
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